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REGULARITY FOR SUBELLIPTIC PDE 
THROUGH UNIFORM ESTIMATES IN MULTI-SCALE GEOMETRIES 


LUCA CAPOGNA AND GIOVANNA CITTI 

Abstract. We aim at reviewing and extending a number of recent results addressing stability 
of certain geometric and analytic estimates in the Riemannian approximation of subRiemannian 
structures. In particular we extend the recent work of the the authors with Rea lfT9l and Manfredini 
ini concerning stability of doubling properties, Poincare’ inequalities, Gaussian estimates on heat 
kernels and Schauder estimates from the Carnot group setting to the general case of Hormander 
vector fields. 
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1. Introduction 


A subRiemannian manifold as a triplet (M, A, ^fo) where M is a connected, smooth manifold 
of dimension n E N, A denotes a subbundle of TM bracket-generating T M, and go is a positive 
definite smooth, bilinear form on A, see for instance [[^ . Similarly to the Riemannian setting, 
one endows (M, A, go) with a metric space structure by defining the Carnot-Caratheodory (CC) 
control distance: For any pair x,y E M set 

do{x, y) = inf{(5 > 0 such that there exists a curve 7 E C°°([0,1]; M) with endpoints x, y 


such that 7 G A( 7 ) and 


Curves whose velocity vector lies in A are called horizontal, their length is defined in an obvious 
way. Subriemannian metrics can be defined, by prescribing a smooth distributions of vector fields 
X = {Xi,, Xm) in K^, orthonormal with respect to go, and satisfying the Hormander finite 
rank condition 


rank Lie{Xi,... ,Xm){x) = n, \/x E 


( 1 . 1 ) 


When attempting to extend known Riemannian results to the subRiemannian setting one nat¬ 
urally is led to approximating the sub-Riemannian metric (and the associated distance function 
do{',')) with a one-parameter family of degenerating Riemannian metric (associated to distance 
functions d^{-, •)), which converge in the Gromov-Hausdorff sense as e —?■ 0 to the original one. 
This approximation is described in detail in from the point of view of the distance functions in Sec¬ 
tion and from the point of view of the Riemannian setting in Definition l3.61 The approximating 
distance functions d^ can be defined in terms of an extended generating frame of smooth vector 
fields Xf, ...,Xp, with p > n and X^ = Xi for i = 1, ...,m, that converges/collapses uniformly 
on compact sets to the original familyXi, ...,Xm as e —)■ 0 . This frame includes all the higher 
order commutators needed to bracket generate the tangent space. When coupled with uniform esti¬ 
mates, this method provides a strategy to extend known Riemannian results to the subRiemannian 
setting. Such approximations have been widely used since the mid-80’s in a variety of contexts. 
As example we recall the work of Debiard Ojl . Koranyi ||5^ [5711 . Ge [|46ll . Rumin fTTl as well 
as the references in Il671l and [[68ll . More recently this technique has been used in the study of 
minimal surfaces and mean curvature flow in the Heisenberg group Starting from the existence 
theorem of Pauls ifTTI . and Cheng, Hwang and Yang lf25]l . to the regularity results by Manfredini 
and the authors ifTSll . [fTbll . Our work is largely inspired to the results of Manfredini and one of us 
[?] where the Nagel, Stein and Wainger estimates for the fundamental solution of subLaplacians 
have been extended to the Riemannian approximants uniformly as e ^ 0. In the following we 
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list in more detail the nature of the stability estimates we investigate. Given a Riemannian mani¬ 
fold (M”, g), with a Riemannian smooth volume form expressed in loeal eoordinates {xi,Xn) 
as d vol = y/gdxi...dxn, one ean eonsider the eorresponding heat operator aeting on funetions 
u : M ^ R, 

1 

LgU = dtu -— ^ di{y/gg^^dju). 

yd ij=i 

The study of sueh operators is elosely related to eertain geometrie and analytie estimates, namely: 
For K CC M and tq > 0 there exists positive eonstants Cd,Cp, .. below depending on K,rQ,g 
sueh that for all t G iT and 0 < r < tq, one has 

• (Doubling property) 

(1.2) vol{B{x,r)) > CDVol{B{x,2r)); 

• (Poineare inequality) \u - UB(x,r)\dvol < Cpr |VgM|ciuo/; 

• (Gaussian estimates) If hg denotes the heat kernel of Lg, x,y E M and f > 0 one has 

(1.3) Cg\vol{B{x, \/t)))~^/^exp(Ag ^^^’^^ ) 

< \h{x,y,t) < Cg{vol{B{x,Vi))) exp(.Bg) 

and if appropriate eurvature eonditions hold 

(1.4) \dldi^...di^h{x,y,t,s) < Cs,kJ~''~Hvol{B{x,t - ); 

L S 

• (Parabolie Harnaek inequality) If LgU = 0 in Q = M x (0, T) and m > 0 then 

(1.5) sup u < Cg inf u. 

B{x,r) X (t-r-2 ,t-r^ /2) B{x,r) X (t+r^/2,t+r^) 

The eonneetions between sueh estimates was made evident in the work of Saloff-Coste fTSll and 
Grigoryan HTll . who independently established the equivalenee 

(Poincare)+(Doubling) <=> Gaussian estimates (11.31) <=> Parabolic Harnaek inequality (11.51) . 
See also related works by Biroli and Moseo |I71, and Sturm llSOll . 

This paper aims at describing the behavior of such estimates along a sequence of metrics g^, 
that collapse to a subRiemannian structure as e ^ 0. ITe will prove that the estimates are stable 
as e -E Q and explore some of the consequences of this stability. Although, thanks to the work 
of Jerison [[5^ . Nagel, Stein and Wainger fTOll and Jerison and Sanehez-Calle [[54ll . the Poineare 
inequality, the doubling property and the Gaussian bounds are well known for subRiemannian 
struetures, it is not immediate that they eontinue to hold uniformly in the approximation as e —)■ 0. 
For one thing, the Riemannian eurvature tensor is unbounded as e —)■ 0, thus preventing the use 
of Li-Yau’s estimates. Moreover, as e —)■ 0 the Hausdorff dimension of the metrie spaces (M, dfj, 
where d^ denotes the distance function associated to g^, typical does not remain constant and in fact 
increases at e = 0 to the homogeneous dimension associated to the subRimannian structure. The 
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term multiscale from the title reflects the fact that the blow up of the metric as e —?■ 0 is Riemannian 
at scales less than e and subRiemannian at larger scales. 

To illustrate our work we introduce a prototype for the class of spaces we investigate, we con¬ 
sider the manifold M = x with coordinates (xi, X2,9). The horizontal distribution is given 
by 

A = span{Xi, X 2 }, with Xi = cos + sin 9dx2, and X 2 = de- 

The subRiemannian metric go is defined so that Xi and X 2 form a orthonormal basis. This is the 
group of Euclidean isometries defined below in example [ATI For each e > 0 we also consider the 
Riemannian metric on M uniquely defined by the requirement that Xi, X 2 , eX^ is an orthonor¬ 
mal basis, with X 3 = — sin 98^^ + cos 98^2 ■ Denote by the corresponding Riemannian distance, 
by X* the adjoint of Xi with respect to Lebesgue measure and by Tg the fundamental solution of 
the Laplace-Beltrami operator Since is uniformly elliptic, then there exists 

Ce, > 0 such that for d^{x, y) < the fundamental solution will satisfy 

C~^d^{x,y)~^ < r,(x,|/) < Ced^{x,y)~^. 


As e —)■ 0 this estimate will degenerate in the following way: ^ ^ 00 and for e 

will eventually have 


To{x,y) ^ do{x,y) 

As a result of the work in llTOl one has that for each e > 0 there exists > 0 such that 

4ix,y) /-p/ dl{x,y) 


C 


-1 


\Be{x,d{x,y))\ 


< r,(x,|/) < C- 


B,{x,d{x,y))\' 


0 one 


The main result of [?] was to provide stable bounds for the fundamental solution by proving that 
one can choose independent of e as e ^ 0. In this paper we extend such stable bounds to the 
degenerate parabolic setting and to the more general subRiemannian setting. 


Since our results will be local in nature, unless explicitly stated we will always assume that 
M = M" and use as volume the Lebesgue measure. The first result we present is due to Rea and 
the authors [|T^ and concerns stability of the doubling property. 

Theorem 1.1. For every cq > 0, and iT CC M” there exist constants i?, C > 0 depending on K, cq 
and on the subRiemannian structure, such that for every e G [0, eo], x E K and 0 < r < R, 

\B^{x,2r)\ < C\B^{x,r)\. 

Here we have denoted by B^ the balls related to the d^ distance function. 

We present here a rather detailed proof of this result, amending some minor gaps in the exposi¬ 
tion in mu. If the subRiemannian structure is equiregular, as an original contribution of this paper, 
in Theorem 13.9! we also present a quantitative version of this result, by introducing an explicit 
quasi-norms equivalent to d^. These families of quasi-norms play a role analogue to the one played 
by the Koranyi Gauge quasi-norm (12.5!) in the Heisenberg group. We also sketch the proof of the 
stability of Jerison’s Poincare inequality from lfT9ll . 
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Theorem 1.2. Let K CC M" and eo > 0. The vector fields (X|)j=i...p satisfy the Poincare in¬ 
equality 


' B^{x,R) 


\U - UB,{x,r)\dx < Cf 


' (a:,2r) 


\Vu\dx 


with a constant Cp depending on K, eo and the subRiemannian structure, but independent of e. 
Here we have denoted by the gradient ofu along the frame Xf X^. 


Our next results eoneems the stability, as e —)■ 0, of the Gaussian estimates for the heat kernels 
assoeiated to the family of second order, sub-elliptic differential equations in non divergence form 

p 

L^^au = dtu-^^ aljX-XjU = 0 , 

Li=i 

in a cylinder Q = O x (0, T). Here {alj}i,j=i,...,p is a constant real matrix such that 

( 1 . 6 ) E «? < Z S 2 A E (I 

2 = 1 = l ^=1 

for all ^ G uniformly in e > 0 and 


(1.7) 


for all ^ G MX and e > 0. 


i=l i,j=i 


2 = 1 


Theorem 1.3. Let K GG M”, A > 0 and Cq > 0. The fundamental solution of the operator 
L^^a, Is a kernel with exponential decay of order 2, uniform with respect to e & [0, eo] and for any 
coefficients matrix A satisfying the bounds above for the fixed A > 0. /n particular, the following 
estimates hold: 


For every K GG Q there exists a constant Ga > 0 depending on A but independent of 
e G [0, Co], and of the matrix A such that for each e G [0, eo], x,y E K and t > 0 one has 


( 1 . 8 ) 


'-'A 


,-Ga 


de{x,yy 


de{x,yr 


< Pe,A^{x,y,t) < Ca- 


g Ca* 

Bfix, ^/t)\ 


\Bfix,^/t)\ 

For s G N and k—tuple (ii,... ,4) G {!,... ,m}^ there exists a constant Cg^k > 0 
depending only on k, s, Xi ,..., A such that 


(1.9) 


\{dlXi^ - ■ ■ Xi^P,^A-){x,yfi)\ < Cg^k 


-2a-k _ME2yr 
t 2 e 

\Bfix, Vi)\ 


for all x,y E K and t > 0. 

• For any Ai, A 2 G M^, s G N and k—tuple {ii,... fik) £ there exists Cg^k > 0 

depending only on k, s, Xi,..., X^, A such that 


WtX,,---x,, 


Pe,A,){x,y,t) - d^Xi^ - ■ ■ Xi^P,^A 2 ){x,y,t)\ < 


( 1 . 10 ) 
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< I l^i — A2\\Cs^k 


h-2s-k 

t 2 e 


\B^{x,Vi)\ 


where 


2 ■= V” a^. 


Moreover, ifVA denotes the fundamental solution of the operator = YlTj=i ihen 

one has 


( 1 - 11 ) ^ • ■■Xi.dlVAo 

e —)■ 0 uniformly on compact sets and in a dominated way on subcompacts ofQ. 

One of our main result in this paper is the extension to the Hormander veetor fields setting of the 
Carnot groups Sehauder estimates established in previous work with Manfredini in [[I]|. To prove 
sueh extension we eombine the Gaussian bounds above with a refined version of Rothschild and 
Stein [TT^ freezing and lifting scheme, adapted to the multi-scale setting, to establish Sehauder type 
estimates which are uniform in e G [0 , Cq], for the family of second order, sub-elliptic differential 
equations in non divergence form 


L^^au ^dtu-'^ alj{x, t)X^X^jU = 0 , 

in a cylinder Q = x (0,T). Our standing assumption is that the coefficients of the operator 
satisfy (11.61) . and (11.71) for some fixed A > 0. 


Theorem 1.4. Let a G (0,1), / G C°°{Q) and w be a smooth solution of L^^aw = f on Q. Let 
K be a compact sets such that K CC Q, set 26 = d^i^K, dpQ) and denote by Kg the 5—tubular 
neighborhood of K. Assume that there exists a constant (7 > 0 such that 


for some value k eN and for every e G [0, eo]. There exists a constant (7i > 0 depending on a, C, 
eo, 6, and the constants in Proposition \5.2\ but independent of e, such that 


re L^fc+2,a/T^\ 6^1 


cfZiKs) 


T IP 


c:v’ {Ks) • 


Here we have set 




\u{x,t) -u{xo,to)\ 


+ sup |m| 


{x,t)jtixo,to) d^{{x,t),{xo,to)) Q 


and ifk>lwe have let u G C^xiQ) if for all i = 1,... ,m, one has Xt G 


k—l,a ( 


Analogous estimates in the spaces, for operators independent of e are well known (see for 
instance fT^ for the constant coefficient case and ifTOll for the Carnot group setting). Our result 
yield a stable version, as e —?■ 0, of such estimates, which is valid for any family of Hormander 
vector fields. 
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Theorem 1.5. Let a G (0,1), / G C°°{Q) and w be a smooth solution of L^^aw = f on Q. Let 
K be a compact sets such that K CC Q, set 25 = do{K, dpQ) and denote by Ks the 5—tubular 
neighborhood of K. Assume that there exists a constant C > 0 such that 

W^ijWcf^iKs) - 

for some value k E N and for every e G [0, eo]. For any p > 1, there exists a constant Ci > 0 
depending on p, a, C, eo, 5, and the constants in Proposition \5.2\ but independent of e, such that 

||u'||^fc+ 2 .p(^) < Cl {\\f\\wf^{Ks) + ■ 

Here we have set 

k 

\Mwf^ 

i=l \I\=i 

2 . Definitions and preliminary results 

Let X = (Xi, denote a colleetion of smooth veetor fields defined in an open subset 

D C M" satisfying Hormander’s finite rank condition ( 11 . 11 ) . that is there exists an integer s such 
that the set of all vector fields, along with their commutators up to order s spans ML for every point 
in VL, 

( 2 . 1 ) rank Lie{Xi,... ,Xm){x) = n, for all Zt G D. 

Example 2.1. The standard example for such families is the Heisenberg group This is a Lie 
group whose underlying manifold is and is endowed with a group law {xi,X2, xf){iii, 1/2,2/3) = 
{xi + 2/1, a;2 + 2/2,2:3 + 2/3 ~ (2^22/1 ~ Xinf)). With respect to such law one has that the vector 
fields Xi = — X 2 dx^ and X 2 = 3^2 + Xid^^ are left-invariant. Together with their commutator 

[Xi,Xf\ = 28x2, yiM a basis ofM^. A second example is given by the classical group of 
rigid motions of the plane, also known as the roto-translation group IZT. This is a Lie group with 
underlying manifoldMf x and a group law (xi, X2, 2/2, ^2) = {xi+yi cos 6 *—2/2 sin 6 , X 2 + 

2/1 sin 0 + 2/2 cos 0 , 01 + 02). 

Following Nagel, Stein and Wainger, [| 70 l page 104 ] we define 

(2.2) xW = {Xi,...,X^}, X( 2 ) = {[Xi,X2],...,[X^_i,X^]}, etc.... 

letting X'^^^ denote the set of all commutators of order k = 1 , ...,r. Indicate by Yi, an 
enumeration of the components of such that Yi = Xi for every i < m. If 

Yk G we say that Yk has a formal degree d{Yk) = d{k) = i. The collection of vector fields 
{Yi,..., Yp} spans at every point. 

Example 2.2. If we consider the Heisenberg group vector fields Xi = — T 2 i 9 x 3 tind X 2 = 

8 x 2 + Xi 8 x 2 with {xi,X 2 ,xfj G then := {Xi,X 2 } and X^^'> = { 203 . 3 }. Ifwe instead 
consider the vectors arising from the group of roto-translations one has Xi = cos 68 x^ + sin 68 x 2 
and X2 = 8 e with {xi,X 2 , 9 ) G x and X ^8 = {Xi,X2} and X^'^^ = {sin0O3;j — cos 003,2}. 
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Example 2.3. Note that the sets may have non-trivial intersection. For instance, consider the 
vector fields 

Xi = cosOdxi + sm.6dx2] X2 = de] X3 = and X4 = 

in {xi,X2, X3, X4, 0 ) G X S^. In this case r = 3 and 

= {Xi,X2,X3,X4}; X(2) = -cos^a,,, 2xsd,,y, and X^^^ = {±Xu2d,,} 

with Yi = Xi,...,Y 4 = X4,Y5 = smOd:^^ - cos 6 * 9 ^ 2,|/6 = 2xzd^^,Yj = Xi,Y^ = -Xi, and 
Yiq = 2dx^. 

2.1. Carnot-Caratheodory distance. For each x,y E Vl and <5 > 0 denote by r(( 5 ) the space of 
all absolutely continuous curves 7 : [ 0 , 1 ] —)■ R”, joining x to y (i.e., 7 ( 0 ) = x and 7 ( 1 ) = y) 
which are tangent a.e. to the horizontal distribution span{Xi, ..., X^}, and such that if we write 

m 

7'(f) = 

i=l 

then^™^ \ai{t)\ < 6 a.e. t G [0,1]. The Camot-Caratheodory distance between x and y is defined 
to be 

(2.3) do(x,y) := inf 6. 

In fTOl . the authors introduce several other distances that eventually are proved to be equivalent to 
do{x, y). The equivalence itself yields new insight into the Carnot-Caratheodory distance. Because 
of this, we will remind the reader of one of these distances. For each x,y E VL and 5 > 0 denote 
by r(( 5 ) the space of all absolutely continuous curves 7 : [ 0 , 1 ] —)■ R", joining x toy and such that 
if one writes 

p 

i=l 

then \ fii{t)\ < One then sets 

d{x,y) := inf 6. 

It is fairly straightforward (see fTOl Proposition 1.1] to see that 

Proposition 2.4. The function d is a distance function in VL and for any K <Z<Z FL there exists 
C = C(Xi,..., Xm, K) > 0 such that 

C-^\x-y\ < d{x,y) < C'|x-2/r""*‘^('). 

It is far less trivial to prove the following (see ITTOI Theorem 4]) 

Theorem 2.5. The distance functions do cind d are equivalent. 
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2.2. The approximating distances. There are several possibile definitions for Riemannian dis- 
tanee funetions whieh approximate a Camot-Caratheodory metrie in the Gromov-Hausdroff sense. 

Definition 2.6. Let {Yi,..., 1^} be a generating family of veetor fields eonstrueted as in (12.21) from 
a family of Hormander veetor fields Xi, For every e > 0 denote by d^{-, •) the Camot- 

Caratheodory metrie assoeiated to the family of veetor fields {Xf, ...,Xp), defined as 

{ Yi if i < m, 

gdW-iy, if rn + 1 < i < p, 

Yi_p+rn i{p+l<i<2p-m 

We will also define an extension of the degree funetion, setting d^{i) = 1 for all i < p, and 
df:{i) = d{i — p + m) if i > p + 1. In order to simplify notations we will denote X = X^, do = d 
and use the same notation for both families of veetor fields (dependent or independent of e). 

Note that for every e G (0, e) the sets {X-} extends the original family of veetor fields (Xj) to 
a new families of veetor fields satisfying assumption (I) on page 107 iTTOll . i.e. there exist smooth 
funetions depending on e, sueh that 

ix-,xa= c‘„xt 

de{l)<de{j)+de{k) 

and 

span W at every point. 

Remark 2.7. Note that the eoeffieients will be unbounded as e —)■ 0. In prineiple this eould be a 
problem as the doubling eonstant in the proof in fTOl depends indireetly from the C"" norm of these 
funetions. In this survey we will deseribe a result, originally proved in [[T9l . showing that this is 
not the ease. 

Remark 2.8. It follows immediately from the definition that for fixed x,y E the funetion de{x, y) 
is deereasing in e and for every e G (0, e), 

do{x,y) > d,{x,y) 

Remark 2.9. Let us eonsider a speeial ease where dim span (Xi,..., X^) is eonstant and the vee¬ 
tor fields Xi, ...,Xp are ehosen to be linearly independent in fl. In this ease we ean eonsider 
two positive defined symmetrie quadratie forms go, and A defined respeetively on the distribution 
H {x) = span (Xi, ...,Xm){x), for t G 12 and on H-^{x). The produet metrie go ® X is then 
a Riemannian metrie on all of Tfl. The form go is ealled a subRiemannian metrie on 12, eorre- 
sponding to H. Next, for every e G (0, e] reeonsider the resealed metrie g^ := go © e~^A and the 
eorresponding Riemannian distanee funetion dg in 12. The latter is bi-Lipsehitz equivalent to the 
distanee d^ defined above. In [|4^ Theorem 1.1] Ge proved that that as metrie spaees, the sequenee 
(12, de) eonverges to (12, do) as e ^ 0 in the sense of Gromov-Hausdorff. In this limit the Hausdorff 
dimension of the spaee degenerates from eoineiding with the topologieal dimension, for e > 0, to 
a value Q > n whieh may ehange from open set to open set. We will go more in detail about 
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this point in the next section. In this sense the limiting approximation scheme we are using can be 
described by the collapsing of a family of Riemannian metric to a subRiemannian metric. See also 
lf 6 ^ Theorem 1.2.1] for yet another related Riemannian approximation scheme. 

Remark 2.10. From different perspectives, note that the subLaplacian associated to the family 
XI, ..., i.e. Cu = is elliptic operator for all e > 0, degenerating to a subelliptic 

operator for e = 0 . 

2.3. A special case: The Heisenberg group EI^ In this section we describe the behavior of the 
distance (and of the corresponding metric balls r) as e —)■ 0 , by looking at the special case 
of the Heisenber group. In this setting we will also provide an elementary argument showing that 
the doubling property holds uniformly as e —?■ 0 . 

Consider the vector fields from Example [2l1 Xi = — X 2 dx ^, X 2 = + Xidx,^ and ^3 = dx^ 

with {xi,X2,x^) G M^. The Camot-Caratheodory distance do associated to the subRiemannian 
metric defined by the orthonormal frame Xi , X 2 is equivalent to a more explicitly defined pseudo¬ 
distance function, that we call gauge distance, defined as 

(2.5) 1x1“^ = (xj + xlf + xl, and p{x,y) = \y~^x\, 

where y~^ = {-yi, -y 2 , -ys) and y~^x = (xi - yi,X 2 - y 2 ,X 3 - ys - {yiX 2 - xiy 2 )) is the 
Heisenberg group multiplication. 

Lemma 2 . 11 . For each x G 

(2.6) < do(a^,0) < A\x\, 
for some constant A > 0 . 

Proof. Observe that the 1-parameter family of diffeomorpthisms 

(xi,X2,X3) (5A(a;i,a;2,a;3) := (Axi, Axa, A^Xg) 

satisfies |dA(a;)| = A|x|, and d5\Xi = XXi o ^a fori = 1,2. Consequently dQ{5\{x),5\{y)) = 
Ado(x, y), and (5A(i?(0,1)) = 5(0, A). Since the unit ball 5(0,1) is a bounded open neighborhood 
of the origin, it will contain a set of the form |x| < A~^ and will be contained in a set of the form 
|x| < A. By applying 6 x we then have that for any 5 > 0, 

{x G M^||x| < A~^R} C 5(0,5) C {x G M^||x| < AR} 

concluding the proof of (12.61) . □ 

Remark 2.12. Since the Heisenberg group is a Lie group, then it is natural to use a left-invariant 
volume form to measure the size of sets, namely the Haar measure. It is not difficult to see OOll 
that the Haar measure coincides with the Lebesgue measure in M^. It follows immediately from 
the previous lemma that the corresponding volume of a ball 5(x, r) is 

(2.7) |5(x,r)| = Cr^ 

As a consequence one can show that the Hausdorff dimension of the metric space (H^, do) is 4. 
The Hausdorff dimension of any horizontal curve (i.e. tangent to the distribution generated by Xi 
and X 2 ) is 1, while the Hausdorff dimension of the vertical z-axis is 2. 
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Next we turn our attention to the balls in the metries and the assoeiated distance functions 
d^. To better describe the approximate shape of such balls we define the pseudo-distance function 
dG,e{x, y) = N^{y~^x) corresponding to the regularized gauge function 

( 2 . 8 ) = xl + xl + min ||x 3 |, 

Our next goal is to show that the Riemannian distance function is well approximated by the 
gauge pseudo-distance dc.e- 

Lemma 2.13. There exists A > 0 independent ofe such that for all x,y 

(2.9) A-^dG,eix, y) < d,{x, y) < AdG,e{x, y). 

The estimate (12.91) yields immediately 

Corollary 2.14. The doubling property holds uniformly in e > 0. 

Remark 2.15. Before proving (12.91) it is useful to examine a specific example: compare two trajec¬ 
tories from the origin 0 = (0, 0, 0) to the point x = (0, 0, Xg). The first is the segment 71 defined 
by s —)■ (0, 0, xgs), for s G [0,1]. The length of this segment in the Riemannian metric given by 
the orthonormal frame Xi, 2 ^ 2 , eXg is 

4(71) = e'VsI- 

We also consider a second trajectory 72 given by the subRiemannian geodesic between the two 
points. In view of (12.61) the length of this curve in the subRiemannian metric g^ defined by the or¬ 
thonormal frame Xi,X 2 is proportional to and coincides with the length in the Riemannian 
metric g"", i.e. 

4(72) = 4(72) ~ 

Since d^ is computed by selecting the shortest path between two points in the g" metric, then if 
V\^\ > e one will have de{x,0) < ^/\^\ ^ N,{ x), whereas at small scales (i.e. for d^^x, 0 ) < e) 
one will have 4(x, 0) < e“^|x 3 |. By left translation invariance of dc.e we have that for any two 
points X = (xi, X 2 , s) and x' = (xi, X 2 , t), 

(2.10) df:{x, x') < C min(e“^|f — s|, \/\t — s|). 

From this simple example one can expect that at large scale (i.e. for points d°(x, 0) > e) the 
Riemannian and the subRiemannian distances are approximately the same d^{x, 0) ~ do{x, 0). 

Proof From the invariance by left translations of both dG,e and d^ it is sufficient to prove that 
d^{x, 0) and N^{x) are equivalent. We begin by establishing the first inequality in (12.91) . i.e. we 
want to show that there exists a positive constant A such that 

A~^N^{x) < 4(0,x). 

Consider a point x = (xi, X 2 , xg) G and three curves 
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A length minimizing curve 7 : [0,1] —)■ for the metric g^, such that 

(ie(0,a:) = 4(7) := [ J «i(^) + ai(f) + e-‘^al{t)dt, 


Jo 

where 7'(t) = Ei=i,3 

• An horizontal curve 71 : [0,1] —)■ with one end-point at the origin (t = 0) and such that 

744 = ai(f)Xi47)+a2(t)X2 4(t). Denote by P = 71(1) and observe that P = (a;i,a;2,P3) 
for some value of such that — p^. 

• A vertical segment 72 : [0,1] —)■ with endpoints P and x, such that 7^4 = ®3(4^3l72(0- 

Note that 


e |a;3-p3|< 




a^{t)dt 


< 


103(4|e dt = 4(72) < 4(7) < d,{x,p). 


Observe that in view of the equivalence (12.61) . 


C ^\jxl + xl < do(^, 0 ) < 4(71) = 4(71) < 4(7), 

for some constant 4 > 0. On the other hand one also has 


e V3 -Psl < de{x,p) < d^{x,0) + d^{0,p) < de{x,0) +4(7i) < 24(0, a;). 

Hence if \ps\ < i|x3| then |x3 — ^3] > i|x3| and consequently 

4 (a;, 0 ) = 4(7) > e'Vs - Psl > min(e”V 3 |, vl^)- 

The latter yields immediately that d^{x, 0) > C~^N^{x), for some value of 4 > 0 independent of 
e > 0. Next we consider the case IpsI > ||x3|. This yields 

min(e"V3 |, a/I^) < ^ min(e“V3|, < VIpsI < 1^1 < Cdo{P, 0) 

< 44(71) = 44(71) < 44(7) = 44 (x, 0 ), 

where |P| is defined as in (12.61) . In summary, so far we have proved the first half of (12.91) . 

To prove the second half of the inequality we consider an horizontal segment Ti joining the 
origin to (^ = (xi, X2, 0). Note that 4(0, Q) = 4(0, Q) = 4(ri) = 4(ri). In view of (12.101) one 
has 

4(0, x) < 4(0, (5) + d^{Q,x) = 4(0, Q) + 4min(e^|x3|, y/jx^) < CN^{x). 

The latter completes the proof of (12.91) . □ 


Remark 2.16. Similar arguments continue to hold more in general, in the setting of Carnot groups. 

As a consequence of Lemma [2.131 one has that for e > 0 the metric space (M^, 4) is locally 
bi-Lipschitz to the Euclidean space, and hence its Hausdorff dimension will be 3. As e —)■ 0 the 
non-horizontal directions are penalized causing a sharp phase transition between the regime at 
e > 0 and e = 0. 

The intuition developed through this example hints at the multiple scale aspect of the d^ metrics: 
At scales smaller than e > 0 the local geometry of the metric space (M^, dj is roughly Euclidean; 
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For scales larger than e > 0 it is subRiemannian. This intuition will inform the proofs of the 
stability for the doubling property in the next seetion. 

3 . Stability of the homogenous structure 

The volume of Camot-Caratheodory balls, and its doubling property, has been studied in Nagel, 
Stein and Wainger’s seminal work ITTOII . In this seetion we reeall the main results in this paper and 
show how to modify their proof so that the stability of the doubling eonstant as e —?■ 0 beeomes 
evident. 

3.1. The Nagel-Stein-Wainger estimates. Consider the Carnot-Caratheodory metrie as- 

soeiated to the family of veetor fields {Xf,X^), defined in (12.41) . Denote by Be{x,r) = 
{y\de{x, y) < r} the eorresponding metrie balls. 

For every n— tuple I = {ii,in) E {1, ...,2p — m}”', and for e > e > 0 define the eoeffieient 

A}(t) = det(X^^(T),...,X^jT)). 

For a fixed 0 < e < e and for a fixed eonstant 0 < (3*2,e < 1, ehoose 4 = (iei,..., im) such that 

(3.1) — ^2,6'n2axjlAj(x)lr‘^‘^'^\ 

where the maximum ranges over all n— tuples. Denote the family of remaining indices, so that 
{Xf^. : e Je} U {X-^ ^ E Je} is the eomplete list Xf ,..., When e = 0 we will 

refer to Iq as a choice eorresponding to the n— tuple X ^^^,..., realizing (13.11) . One of the main 

eontributions in Nagel, Stein and Wainger’s seminal work [ITOll . eonsists in the proof that for a v 
and a x fixed, and letting 

Q,(r) = {uEW : \uj\ < 

denote a weighted eube in M”, then the quantity |Aj^(t)| provides an estimates of the Jaeobian of 
the exponential mapping u ^e,v,x{u) defined for u E Q{r) as 

( 3 . 2 ) <I>,,„,^(m) = ea;p(^ ^ 


More preeisely, for e > 0 and fixed one has 
Theorem 3.1. IITOl Theorem 7] 

For every e > 0, and K CG MF there exist > 0 and constants 0 < (7i e, C 2 ,e < 1 such that 
for every x E K and 0 < r < R^, if R is such that (13.11) holds, then 

i) if\Vk\ < is one to one on the box (5e(Ci,er) 

ii) if\vk\ < the Jacobian matrix of^^^y^x satisfies on the cube QfiCi^^^r) 

7 |Al^(x)| < \ J^e,v,x\ < 4|A^^(x)| 


iii) 


^e,vAQe{Cyer)) C Bfix,r) C ^e,v,x{Qe{Cyer/C2,e)) 
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As a corollary one has that the volume of a Carnot-Caratheodory ball centered in x can be 
estimated by the measure of the eorresponding eube and the Jaeobian determinant of 

Corollary 3.2. l llTOl Theorem l]j For every e > 0, and K CG'MF and for > 0 as in Theorem 
I3.il there exist constants C^ei > 0 depending on K, R^, and C 2 e such that for all x E K 
and 0 < r < Rf^ one has 

(3.3) < \B^{x,r)\ < ^4, ^ 

I I 

Estimates (13.31) in turn implies the doubling eondition (11.21) with eonstants depending eventually 
on Cu and C 2 e- 


3.2. Uniform estimates as e —)■ 0. Having already proved the stability of the doubling property in 
the speeial ease of the Heisenberg group, in this seetion we turn to the general ease of Hormander’s 
veetor fields and deseribe in some details results from [[T^ establishing that the eonstants Cu C 2 e 
do not vanish as e —)■ 0. Without loss of generality one may assume that both eonstants are 
non-deereasing in e. In faet, if that is not the ease one may eonsider a new pair of eonstants 
Q,, = inf^g[,,,-] Ci^s, for i = 1, 2. 

Proposition 3.3. For every e e [0, e], the constants R^, and € 2 ,^ in Theorem \3.1\ mav be chosen 
to be independent ofe, depending only on the norm of the vector fields, on the number e, and 

on the compact K . 


Proof. The proof is split in two oases: First we study the range e < r < Rq whieh roughly 
oorresponds to the balls of radius r having a sub-Riemannian shape. In this range we show that 
one oan seleet the eonstants to be approximately Cifl. The seoond ease eonsists in the analysis 
of the range r < e < e. In this regime the balls are roughly of Euolidean shape and we show that 
the eonstants (^'j ^ ean be approximately ohosen to be Ci^i. 

Let us fix e G (0, e], i? = i?o and r < Rq. We oan start by describing the family R defined in 
(13.11) . which maximize A^x). We first note that for every e > 0 and for every i,m + l<i<pwe 
have 


(3.4) 




de{i+p-m) _ d{i)Y. 


In the range 0 < r < e < e one oan assume without loss of generality that the n—tuple satisfying 
the maximality eondition (13.11) will inolude only veetors of the form 

for some n—index R = (f^i, ...,ien), with I < iek R P- In fact, if this were not the ease and 
the n—tuple were to inolude a veetor of the form Xj = Yj_p^rn for some p < j, then we oould 
substitute suoh veetor with (|3.4[) infer that the value of 

the eorresponding term \Xf(x)\r^dtc) would inerease. 

Similarly, in the range 0 < e < r < e one oan assume that the n—tuple satisfying the max¬ 
imality eondition (13.11) will inelude only veetors of the form for some n—index 

h = {Ri, ■■■fieri), with I < iek R ?■ Noto that the eorresponding expression 




.d,(h)-l _ 


= |det(ri^,,...,yi^J(a:)|r 
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would then be one of the terms in the left hand side of (13.11) for e = 0, and thus is maximized by 

Case 1: In view of the argument above, for every e < r < Rq the indices R defined by the 
maximality condition (13.11) can be chosen to coincide with indices of the family R and do not 
depend on e. On the other hand the vector excluded from R will be not only those in Jq but also 
the ones that have been added with a weight factor of a power of e, 

{XI-.ke R} = : to,k e Jo} u ^ : to,k e Io,to,k > m} 

In correspondence with this decomposition of the set of indices we define a splitting in the t>—variables 
in (15.141) as 

V = {v, V, v). 

Consequently for every e < r the function ^ Ru) can be written as 
(3.5) 

= exp'^kXl^Yx) = expY = 

expl^Y RKj+ ^ Y + Y ( t ) = 

*OjS/o iok^Jo *0fcS/0i*>™' iok^Jo,io,k>nT- 

Let us define mappings 

Fl,eA^) = •••, Wm, Um+1 + ..., , 

and 

^2,.(w) = ...,V 2 p-m + • 

In view of (13.51) we can write 


(3.6) 


(^) ®0,T2,e(D),a; (^) ) • 


Note that for any e > 0 and for a fixed v, the mapping u Fi^^^y{u) is invertible and volume 
preserving in all M". Moreover = J$o,F 2 «(«)). In view of (13.61) and of Theo¬ 

rem [3T1 as a function of u, the mapping ^e,v,x{u) is defined, invertible, and satisfies the Jacobian 
estimates in Theorem 13.1 1 (ii) 


1 

4 



< 


J*^0,T2,e(^’),a; (-^1,6,1)(^)) I 


M,,yM\ < RXl{x)\ 


for all u such that Fi ^ „(«) G Qo{Cifir) and for v such that 


FY^)\ = 14 + 


4i| < • • • \uR < \um+i + 

when A; = 1,..., 2 p — m. 
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The completion of the proof of Case 1 rests on the following two claims: 

Claim 1 let e < r < Rq. There exists Cq > 0, independent of e, such that for all v satisfying 
\vk\ < one has < ( 72 , 

Proof of the claim: If we choose Cq < min{(7i,o, < 72 , 0 } and 


\vk\, \vk\, l^fcl < min{(7i,o, ( 72 ,o}^^, \uj\ < (7i,o 


it follows that 


fpd{iQk') rp 

\vk\ < < 72,0 ^ , l^fcl, \vk\ < <7i,o^, l^il < <7i,o 


4 ’ 


So that 


^ (72,0" 


4 ’ 




< ^1,0^—, |m,| < Cl, 


^d{ioj ) 

4 ’ 


completing the proof of the claim. 

Claim 2 Let e < r < Rq and v fixed such that \vk\ < (76r'^7efc) for k = 1,2p — m. One has 

that 

Q.(C5-V) C Fi;,^,(Qo(Ci,or)) C Q.{C,r) 

for some constant (75 > 0 independent of e > 0. 

Proof of the claim: Choose sufficiently large so that 2max{(77^, Cq} < (7i,o and observe 
that if M G Q^{C'^^r) then for /c = 1, ...,m we have \uk\ < (7i,or'^<*'’'“^ = (7i,or’'^<*“’''^ while for 
A: = m + 1, n we have ,,(«)! = \uk + VkC^'‘°^C^\ < max{( 77 ^, ( 76 }r'^<®o'“)(l + e'^7ofe)-i) ^ 
( 7 i or'^(®ofc). This proves the first inclusion in the claim. To establish the second inclusion we choose 
Cq large enough so that 2((7i,o + ( 72 ,e) < Cq and observe that if Fi^^^^{u) G QaiCi^r) then for 
k = m + 1, ..., none has \uk\ < \uk + VkC‘^'^°^"^~^\ + \vk\C^'^°^C^ < 2((7i,o + ( 72 ,e)r'^(*°'“) < CqC^^°*‘\ 
The corresponding estimate for the range A; = 1,..., m is immediate. 

In view of Claims 1 and 2, and of Theorem 13.11 It follows that for e < r and these choices of 
constants (independent of ejilthe function ^e,v,x{u) is invertible on Qo{Ciflr) and i), ii) and iii) are 
satisfied. 


Case 2: As remarked above, in the range 0 < r < e < e one can assume that the n—tuple satisfy¬ 
ing the maximality condition (13.11) will include only vectors of the form ..., e^Cn)-^Yi^^} 

for some tt,— index R = (i^i, •••, Rn)^ with 1 < < p. Note that in view of (13.41) and the maxi¬ 

mality condition (13.11) the corresponding term 

can be rewritten and estimated as follows 

|A',. Wl/A'-) = e‘‘<'-)-"r"| det(y.,.,y..J(a:)|. 


*i?o in place of R^, Cq in place of Ci,e and Cq in place of C2,e 
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It is then clear that the maximizing n—tuple 4 in (13.11) will be identified by the lowest degree d{I^) 
among all n—tuples corresponding to non-vanishing determinants det(yi^^, Yi^^) in a neighbor¬ 
hood of the point x. Since this choice does not depend on e > r, then one has that 4 = I^. In other 
words, if we denote 


r — ... \ 


then the maximality condition (13.11) in the range 0 < r < e < e can be satisfied independently 
from e by selecting the family of vector fields: 




The complementary family becomes 


(3.7) e JJ 

= : Zo,fc G Ji, with < p} U : i^^k G J^, with > p} 

If we denote A^, and these three sets, and split the n—variable from (15.141) as n = {v, v), then it 
is clear that 

and in this case the values of and d^ are the same on the corresponding indices. Analogously 
Y G iff y G and the degrees are the same. 

For every e > r the map then can be written as 


^e,v,xH = exp ( = e^P ( “7^4 + Y1 


= exp 


E 




fcGt/e 


7e? 


iek 


[X] 


u 


xl 


7 


E 


c:'7(*0.fc)-l 


Vk- 


iek^Je and iej<p 

This function is defined and invertible for 




A + E 


iik&Je and Uj>p 


\X) 






U 






Recall that with the present choice of r < e < e, we have = (3'i,er'^'P'^7) = Ci 

If we set 

I, |wfc| < 


\Ui 




and argue similarly to Case 1, then the function ^e,v,x will satisfy conditions i), ii), and iii) on 
Q{Ci^ir) and hence on QiCi^f^r), with constants independent of e. □ 
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3.3. Equiregular subRiemannian structures and equivalent pseudo-distances. The intrinsic 
definition, based on a minimizing choice, of the Camot-Caratheodory metric is not convenient 
when one needs to produce quantitative estimates, as we will do in the following sections. It 
is then advantageous to use equivalent pseudo-distances which are explicitly defined in terms of 
certain system of coordinates. In the last section we have already encountered two special cases, 
i.e. the norms | ■ | defined in (12.51) and its Riemannian approximation (12.81) . In this section we 
extend this construction to a all equi-regular subRiemannian structures. For C M" consider the 
subRiemannian manifold {Q,, A, g) and iteratively set := A, and A*+^ = A* + [A*, A] for 
i G N. The bracket generating condition is expressed by saying that there exists an integer s G N 
such that Ap = M" for all p G M. 

Definition 3.4. A subRiemannian manifold (fl, A, g) is equiregular if, for all i G M, the dimension 
of Ap is constant in p G 12. The homogenous dimension 

s-l 

(3.8) Q = ^[dim(A;+^) - dim(A*)], 

i=l 

coincides with the Hausdorff dimension with respect to the Carnot-Caratheodory distance. 

This class is generic as any subRiemannian manifold has a dense open subset on which the 
restriction of the subRiemannian metric is equiregular. 

Example 3.5. Systems of free vector fields, as defined in Deiinition \5.4\ yield a distribution A that 
supports an equiregular subRiemannian structure for any choice of the horizontal metric g. 

Next we assume we have a equiregular subRiemannian manifold (12, A, g) and consider an or¬ 
thonormal horizontal basis Xi, of A. Following the process in (12.21) one can construct a 

frame Yi,..., for M" where Yi,..., Ym is the original horizontal frame and Y^+i,..., Y„ are com¬ 
mutators such that (Yi,..., Ym^)\p spans A^, for A; = 1,..., s. The degree d{i) of Yj is the order 
of commutators needed to generate Yj out of the horizontal span, i.e. d{i) = fc if Yj G A^ but 
Yj ^ Ap~^. In particular one has d{i) = 1 for i = 1,..., m. The equiregularity hypothesis allows 
one to choose Yi,..., Y„ linearly independent. Next we extend p to a Riemannian metric gi on all 
of T12 by imposing that Yi,..., Y„ is an orthonormal basis. 

Definition 3.6. For any e G (0, e] we define the Riemannian metric g^ by setting that Yj, 

i = 1,..., n is an orthonormal frame. Denote by dfix, y) the corresponding Riemannian distance 
function. 

Remark 3.7. Repeating the proof of IITOl Theorem 4] one immediately sees that d^ as defined here 
is comparable to the distance defined in Section with equivalence constants independent of 

e > 0. 

We define canonical coordinates around a point Xq G 12 as follows. Since Yi,..., Y„ is a generat¬ 
ing frame for T12 then for any point x in a neighborhood ce of Xq one has that there exists a unique 
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n—tuple (xi, Xn) such that 

n 

(3.9) exp(^a:iy))(a;o) 


= X. 


i=l 


We will set x = (xi, ...,Xn) and use this n—tuple as local coordinates in cj. 

Definition 3.8. For every x = (xi,...,Xn) G ce we define a pseudo-distance dG,e(x,xo) := 

JV^(xi, ...,Xn) with 


(3.10) N^{xi,...,Xn) : = 

For e = 0 we set 

No{xi, ...,Xn) : = 


xf + min (e 

\| i=i 


I ? \^i 


\l/d{i) 


)■ 


2=m+l 


\ 


2=1 2=m+l 


\X, 




Theorem 3.9. For every compact xq E K (Z u there exists C = C{K, A, g, u) > 0, independent 
o/e G (0, e], such that 

C~^dG,e{x,XQ) < de{x,Xo) < CdG,e{x,Xo) 

for all X E K. 


Remark 3.10. Note that for e = 0 the equivalence is a direct consequence of the Ball-Box theorem 
proved by Nagel, Stein and Wainger |[70]l or Mitchell ll65l Lemma 3.4]. This observation replaces 
the estimates (12.61) from the Heisenberg group setting. 

The proof of Theorem follows as a corollary of the following 

Proposition 3.11. In the hypothesis ofTheorem \3.9\ one has that there exists R = R{K, A, g, u) > 
0,C = C{K, A, g, u) > 0, independent ofe E (0, e] , such that for all x E K and r E (0, R), 

BG,e{xQ,C~^r) C 5^(xo,r) C BG,^{xQ,Cr), 

where 





M" such that 


max 

2 = 1,...,S 


min 



Proof The proof follows closely the arguments in the previous section and is based on the results 
in irTOll . In view of the equiregularity hypothesis note that Yi,..., are linearly independent and the 
construction in (12.41) yields the distribution Xf X^n-m Recall from (15.141) . Proposition 

13.31 and Theorem |3.II that if R, are chosen as in (13.11) and for any v = (ui,..., Vn-m) such that 
\vk\ < one has 

(3.11) BRxQ.r) ^ ^e,v,xo{Qe{r)), 
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with constants independent from e > 0, where Qe = {m G : \uj\ < and 

^e,vAu) = exp(^ UjXl .+ Y ^kXl^^Yx). 

The n— tuple 4 contains n indexes related either to the horizontal vector fields XI, ■■■,X^ or 
to the eommutators X^_^_^,X^. The latter may eonsist of weighted versions ...,X^ 

or unweighted versions In either ease the same vector will appear both in the 

weighted and in the unweighted version (either among the indexes or in the eomplement J^). 
Comparing the representation with the x—eoordinates representation (13.91) one has 

n 

exp(^Xiri)(xo) = Y ^kXl^^{xo), 

i —1 

and we let for eaeh k = 1, ...,n 

gd(fc)-i^.^ + Vj^ ifik<n 
Ui^ + if 4 > n ■ 

From the latter we obtain that for all fc = 1, n 

\xk\ < + 

If x G B^{xo,r) then \uij\, \vj^\ < Consequently, 


^d{k)-l^ ^d{k) 


1 / aytv) 



< C minfr 


i.i-j 


\ d{k—l)- 


\ d(k)-l 

-) 

r ) 


l/d(/c) 


) < 2Cr. 


This shows that for r > 0 sufficiently small, and for some ehoiee of (7 > 0 independent of e > 0, 
we have B^{xq, r) C BG,t{xo, Cr). 

To prove the reverse inclusion we consider a point x = exp(^”^^ XiYi}{xo) G BG,e{xo, Cr). 
Select 4 as in (13.11) and set u = 0 to represent x in the basis ..., Xi^ as 


x = exp{^ Y UjXl^^{xo). 

In view of Theorem 13.1[ and (13.1 II) . to prove the proposition it suffiees to show that there exists a 
eonstant (C > 0 independent of e > 0 sueh that for each j = 1,..., n one has \uj\ < Cr‘^Xej)_ 

We distinguish two oases: In the range e > 2r one can argue as in (13.41) to deduce that for each 
j = 1,..., n we may assume without loss of generality that the contribution due to UjX^^ . follows 

from the choice of a weighted veetor, and hence is of the form for some k > m. 

Consequently one has d^{i^^j) = 1 and Xk = 

On the other hand, sinoe e > 2r then one must also have that 
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Consequently one has 

\uj\ = = 

In the range e < 2r we observe that one must have \xk\ < Cr^{k). Arguing as in (13.41) we see 
that without loss of generality, or eaeh j = 1, the eontribution due to UjXf^. follows from 
the ehoiee of a un-weighted veetor, and henee is of the form UjYk for some k > m. Consequently 
one has de{ij) = d{k) > 1 and Xk = Uj, eoneluding the proof. 

□ 


4 . Stability of the Poincare inequality 


In this seetion we will foeus on the Poineare inequality and prove that it holds with a ehoiee of a 
eonstant whieh is stable as e —)■ 0. Our argument rests on results of Laneonelli and Morbidelli [(601 
whose proof, in some respeets, simplifies the method used by Jerison in Il5^ . Using some Jaeobian 
estimates from [l45ll or iHTll we will establish that the assumptions required in the key result ll60l 
Theorem 2.1] are satisfied independently from e > 0. We start by reealling 


Theorem 4.1. Theorem 2.1] Assume that the doubling condition (D) is satisfied and there 
exist a sphere Bfxo,r), a cube C M” and a map E : Bfxo,r) x ^ MA satisfying the 
following conditions: 

i) Bfxo,2r) C E{x,Qfj for every x E Bfxo,r) 

ii) the function u i-a E{x, u) is one to one on the box as a function of the variable u and 
there exists a constant gi > 0 such that 

— \JE{x,0)\ < \JE{x,u)\ < ai\JE{x,0)\ for every u E 

Oil 


Also assume that there exists a positive constant a 2 , and a function 7 : BfxQ, r) x Q^x [ 0 , a 2 r] -E- 
M” satisfying the following conditions 

iii) For every (t, u) E Bfxo, r) x the function t ha 7 (x, m, t) is a subunit path connecting 
X and E{x, u) 

iv) For every {h, t) E Bfxo, r) x the function x ha y[x, u, t) is a one-to-one map and there 

exists a constant ^3 > 0 such that 


inf 

iJe(xo,r)xQe 


det 


dy 

dx 


> 03 


Then there exists a constant Cp depending only on the constants Oi, 02 , 03 cind the doubling con¬ 
stant Cp such that (P) is satisfied. 


We are now ready to prove Theorem 1 1.21 

Proof All one needs to establish is that the assumptions of Theorem 14. II are satisfied unformly in 
e on a metrie ball. Apply Proposition 13.31 and Theorem 13.11 with K = Bfixo,r) and choose the 
constants Ci produced by these results. Set = Qfi^r) and let 

E{x, u) = ^e,o,xiu), defined on K x ^ M". 
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To establish assumption (i) of Theorem 14. II it suffiees to note that by virtue of eondition (hi) in 
Theorem [3T] one has that for x G i?e(xo, r), 


B^{xo,2r) C B^{x,3r) C E{x,Q^). 

Assumption (ii) in Theorem 14.11 is a direet eonsequenee of eondition (ii) in Theorem I3.1[ with 
= 16. Chow’s eonneetivity theorem implies that E{x,u) satishes assumption (hi), with a 
funetion 7 , pieeewise expressed as exponential mappings of veetor helds of e—degree one. Let us 
denote the required veetor helds. With this ehoiee of path, it is known (see for example 

Il45l Lemma 2.2] or ||4T1 pp 99-101]) that x —)■ 7 (x, u, t) is a path, with Jaeobian determinant 


d'y 

det-;—(x, u, t) 
ox 


1 + i>{x,u,t), 


for a suitable funetion 'ip^x, u, t) satisfying 

|'^(x, u,t)\ < cr, on K x x [ 0 , cr]. 

Sinee the eonstant c depends solely on the Lipsehitz eonstant of the veetor helds then it 

ean be ehosen independently of e. As a eonsequenee eondition (iv) is satished and the proof is 
eoneluded. 

□ 


5 . Stability of Heat Kernel Estimates 

5.1. Hormander type parabolic operators in non divergence form. The results in this seetion 
eoneem uniform Gaussian estimates for the heat kernel of eertain degenerate parabolie differential 
equations, and their parabolie regularizations. We will eonsider a eoheetion of smooth veetor helds 
X = (9fi, • • • , Xm) satisfying Hormander’s hnite rank eondition (11.1!) in an open set 12 C We 
will use throughout the seetion the dehnition of degree d{i) relative to the stratiheation ( 12 . 21 ) . 

A seeond order, non-divergenee form, ultra-parabolie operator with eonstant eoefheients Ojj ean 
be expressed as: 

m 

( 5 . 1 ) LA = dt-Y,^^ 3 XiXj, 

where A = (ajj)jj=i is a symmetrie, real-valued, positively dehned m x m matrix satisfying 

m 

(5-2) A-' E ^ E ^ E 

d{i) = l hj = ^ d{i) = l 

for a suitable eonstant A. We will also eall 

(5.3) Mm,A the set of symmetrie m x m real valued matrix, satisfying (15.21) 

If A is the identity matrix then the existenee of a heat kernel for the operator La is a by now 
elassieal result due to Folland Il38]l and Rothsehild and Stein fT^ . Gaussian estimates have been 
provided by Jerison and Sanehez-Calle ll54ll . and by Kusuoka and Strook ll59ll . There is a broad, 
more reeent literature dealing with Gaussian estimates for non divergenee form operators with 
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Holder continuous coefficients Uij. Such estimates have been systematically studied in |l8l, [fTOl . 
[HI, [fT^ where a self-contained proof is provided. 

A natural technique for studying the properties of the operator La is to consider a parabolic 
regularization induced by the vector fields Xf defined in (12.41) . More precisely, we will define the 
operator 

p 

(5.4) = 

where ajj is any p x p positive defined matrix belonging to Mp 2 A and such that 

= forz, j = 1,.. .m. 

We will denote 

(5.5) 

the set of such matrices. Formally, the operator La can be recovered as a limit as e —)■ 0 of operator 
Lg A- Here we are interested in understanding which are the properties of solutions of L^^a which 
are preserved in the limit. 

For e > 0 consider a Riemannian metric defined as in Remark IZ9l such that the vector fields 
XI are orthonormal. The induced distance function is biLipschitz equivalent to the Euclidean 
norm 11^;. Consequently, the operator L^ a has a fundamental solution F^^^, which can be estimated 
as 

(5.6) < c.^-^ 

for some positive constant depending on A, e and Xi,X^. 

Unfortunately the constant blows up as e approaches 0, so the Riemannian estimate (15.71) 
alone does not provide Gaussian bounds of the fundamental solution F^ of the limit operator (15.11) 
as e goes to 0. In [[58ll the elliptic regularization technique has been used to obtain L^ and C°‘ 
regularity of the solutions, which however are far from being optimal. In [|28]|, new estimates 
uniform in e have been provided, in the time independent setting which are optimal with respect to 
the decay of the limit operator. In [[T71l the result has been extended to the parabolic operators, in 
the special case of Carnot groups. 


In order to further extend these estimates, we need to formulate the following definition: 

Definition 5.1. We say that a family of kernels {Pe,A)t>o,A&M^ 2 a’ defined on oo[ has, on 

the compact sets of an open set U, an exponential decay of order 2 + h, uniform with respect to 
a family of distances {d^)e and of matrices A G Mp 2 A (see definition 15.51) and we will denote 
Pe,A E S{2 + h, de, Mp 2 a) if the following three condition hold: 
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(5.7) 


(5.8) 


For every K <Z<Z Vt there exists a eonstant Ca > 0 depending on A but independent of 
e > 0, and of the matrix A G 2 a sueh that for eaeh e > 0, x,y ^ K and t > 0 one has 



tie t 

\B^{x, V^)| 


< Pe,Aix,y,t) < Ca 


h dt(x,y)'^ 

tie 


B^{x, '/t)\ 


For s G M and /c—tuple (ii,..., 4) G there exists a eonstant Cs,k > 0 de¬ 

pending only on /c, s, Xi,A sueh that 




^ Cs,k 


h- 2 s-k 

tie 


B^{x,Vt)\ 


for all x,y E K and t > 0. 

For any Ai, A 2 G Ma, s G N and /c—tuple (A,..., 4) £ there exists Cs,k > 0 

depending only on fc, s, Xi,..., X^, A sueh that 


(5.9) 


\{dlX,, • • ■ X,,P,,aJ(x, y, t) - dtXi, • • ■ X,,P.,aJ(x, yA)\< 


< IIAi-AallC',,^ 


h-2s-k _ de(x,y) 

tie '^A* 

\B^{x,Vi)\ 


where ||A||2 := 


With these notations we will now extend all these previous results to veetor fields whieh only 
satisfy the Hormander eondition, establishing estimates whieh are uniform in the variable e as 
e —0, and in the ehoiee of the matrix A G for the fundamental solutions T^ a of the operators 
L^^a- To be more speeifie, we will prove: 


Proposition 5.2. The fundamental solution T^^a of the operator L^^a, is a kernel with exponential 
decay of order 2, uniform with respect to e > 0 and to A E according to definition i l5.iD . 

Hence it belongs to the set S{2,de, M^a)- Moreover, z/Ta is the fundamental solution of the 
operator La defined in di.il) one has 

(5.10) X/^ • --XldlV^A ^ • ■■Xi.dlVA 

as e ^ 0 uniformly on compact sets and in a dominated way on subcompacts ofVL. 


Our main eontribution is that all the eonstants are independent of e. The proof of this assertion is 
based on a lifting proeedure, whieh allows to express the fundamental solution of the operator La,^ 
in terms of the fundamental solution of a new operator La independent of e. The lifting proeedure 
is eomposed by a first step in whieh we apply the delieate Rothsehild and Stein lifting teehnique 
fy^ . After that, when the veetor fields are free up to a speeifie step, we apply a seeond lifting 
whieh has been introdueed in Il28]l . where the time independent ease was studied, and from ifTTll 
where the Carnot group setting is eonsidered. 

The simplest example of sueh an equation is the Heat equation assoeiated to the Kohn Laplaeian 
in the Heisenberg group, Lt — X^ — X|, where the veetor fields Xi and X 2 have been expressed on 
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coordinates in Example 12. II In order to present our approaeh we will give an outline of the proof 
in this speeial setting. 

Example 5.3. Denote by {xi,X 2 ,x^) points ofM.^, let Xi,X 2 ,X 3 be the vector fields defined in 
Example \2.1\ and let I denote the identity matrix: Consider the parabolic operator 

Uj = -dt + xl + xl + ^xl 

and note that it becomes degenerate parabolic e —)■ 0. Let d^ denote the Carnot-Caratheodory 
distance associated to the distribution Xi, X 2 , 6 X 3 . 

In order to handle such degeneracy we introduce new variables (zi, ^ 2 , -^ 3 ) cind a new set of 
vector fields replicating the same structure of the initial ones, i.e., 

Z\ dzi ~\~ Z 2 dz 3 1 Z 2 Qz 2 9 z 3 ) dz 3 

with (ti, X2, X3, zi, Z 2 , Z 3 ) G X EI^ The next step consists in lifting to an operator 

Le = dt + Xi + X 2 + Z‘^ + Z 2 + (Z 3 + eXs)^, 

defined on x and denote by its fundamental solution. Let d^ denote the Carnot- 
Caratheodory distance generated by Xi,X 2 , Zi, Z 2 , {Z 3 + eXs) and arguing as in (15.221) note 
that d,.{{x, z), {y, z)) > de{x,y)—C q, for some constant Cq independent of e. Consider the change 
of variables on the Lie algebra o/H^ x 


Xj —)■ Xj, Zi —)■ Zi, for i = 1 , 2 , Z 3 + 6 X 3 —)■ Z 3 . 


Note that the Jacobian of such change of variables does not depend on e and that it reduces the 
operator to 

L = dt -\- X^ + X 2 + Z"^ + Z 2 + Z"^ 

whose fundamental solution we denote by f. Note that this operator is parabolic with respect to 
the vector fields Zi and degenerate parabolic with respect to the vector fields Xj. Is is clear that 
the operator L is independent of e, and consequently its fundamental solution f satisfies standard 
Gaussian estimates with constants independent of e 


r(x,t) < Cx 


d(x,0)^ 
g Cxt 




where d denotes the Carnot-Caratheodory distance in x generated by the distribution of 
vector fields Xi, X 2 , Zi, Z 2 , Z 3 . Changing back to the original variable we see that also f ^ sat¬ 
isfies analogous estimates with the same constants, with the distance d replaced by the distance 
de naturally associated to the operator Linally, integrating with respect to the added variable 
[zi, Z2, Z3), we obtain an uniform bound for the fundamental solution of the operator in terms 

of the distance d,.. 
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5.2. The Rothschild-Stein freezing and lifting theorems. Let us first recall a local lifting pro¬ 
cedure introduced by Rothschild and Stein in ITT^ which, starting from a family of 

Hormander type vector fields of step s in a neighborhood of leads to the construction of a new 
family of vector fields which are free, and of Hormander type with the same step s, in a neighbor¬ 
hood of a larger space. The projection of the new free vector fields on yields the original vector 
fields, and that is why they are called liftings. 

Let us start with some definitions: 

Definition 5.4. Denote by nm,s the dimension (as a vector space) of the free nilpotent Lie algebra 
with m generators and step s. Let Xi,, Xm be a set of smooth vector fields defined in an open 
neighborhood of a point xq e M”, and let 

= span{X^^\--- ,Xh)}, 

where the sets X^ are as defined in (12.21) . We shall say that Xi,..., X^, axe free up to step s if for 
any 1 < r < s we have nm,s = dim{y^^'>). 

If a point Xq G i?" is fixed, the lifting procedure of Rothschild-Stein locally introduces new 
variables 2 and new vector fields {Zfj expressed in terms of the new variables such that in a neigh¬ 
borhood U of Xq the vector fields Xi = {Xi + are free at step s. More precisely, one has 

fy^ Theorem 4] 

Theorem 5.5. Let Xi,..., Xm be a system of smooth vector fields, satisfying (li.il) in an open set 
U C For any x E U there exists a connected open neighborhood of the origin V C and 

smooth functions Xij{x, z), with a; G i?" and z = {zn+i,... ,zf) G V, defined in a neighborhood 
U ofx = (x, 0) G f/ X L C R^, such that the vector fields Xi,..., X^ given by 

V 

Xi=XiF Zi, Zi= ^ \ij{x, z)d^. 

j=n+l 

are free up to step r at every point in U. 

Remark 5.6. In the literature the lifting procedure described above is often coupled with another 
key result introduced in fi^ . a nilpotent approximation which is akin to the classical freezing 
technique for elliptic operators. Let us explicitly note that in section 5.3 we need only to apply the 
lifting theorem mentioned above, and not the freezing procedure. In particular, in the example of 
the so called Grushin vector fields 

^3 = 5x1 and X 4 = xida:^ 

they would need to be lifted through this procedure to the Heisenberg group structure 

X3 = and X4 = -f xA^. 

On the other hand the vector fields 

Xi = cos -f sin 68 x 2 and X2 = dg 
will be unchanged by the lifting process, since they are already free up to step 2. 
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Later on, In section 5.4 we will apply Rothschild and Stein’s freezing theorem to a family of 
vector fields Xi,..., free up to step r. This will allow to approximate a given family of vector 
fields with homogeneous ones. Note that in this case the function $ in (15.141) is independent of v 
and its expression reduces to: 

(5.11) = exp(y^^UiXj){x). 

i 

The pertinent theorem from is the following, 


Theorem 5.7. Let Xi,... ,Xm be a family of vector fields are free up to rank r at every point. 
Then for every x there exists a neighborhood V of x and a neighborhood U of the identity in Gm,r> 
such that: 

(a) the map '■ U V is a dijfeomorphism onto its image. We will call Qx ds inverse map 

(b) we have 

(5.12) dQx{Xi) = Yi + Ri, i = l,...,m 

where Ri is a vector field of local degree less or equal than zero, depending smoothly on x. 
Hence the operator Ri will represented in the form: 

Ri = '^ai{u)Xi, 
jh 

where a is an homogeneous polynomial of degree d{Xi) — 1. 


5.3. A lifting procedure uniform in e. So far we have started with a set of Hdrmander vector 
fields Xi, ..., Xm in C and we have lifted them through Theorem [53] to a set Xi, ..., X^ of 
Hormander vector fields that are free up to a step s in a neighborhood kl C ML. Next, we perform a 
second lifting inspired by the work in ifTTll . We will consider the augmented space x defined 
in terms of u new coordinates 2 = {zi, ...,2^). Set .2 = (5, z) and denote points of x W by 
X = (x, 5, z) = (x, z). Denote by Zi, ..., a family of vector fields free up to step s. Xi, ..., Xm, 
i.e. a family of vector fields free of step s in the variables z, and let 

Zm+l ) ■ ’ ’ 


denote the complete set of their commutators, as we did in (12.21) . Note that the subRiemannian 
structure generated by Zi,..., Z^ coincides with the structure generated by the family Xj, but are 
defined in terms of new variables z. 

For every e E [0,1) consider a sub-Riemannian structure determined by the choice of horizontal 
vector fields given by 


(5.13) (X^, • • • X^_,_J,) — (Xi, . . . , Xm, Zi, . . . , Zm, X^_,_^ + Zm+1, . . . , X^ + Zj,). 

Since the space is free up to step r the function $ in ( 15 . 141 ) is independent of v and its expression 
reduces to: 

( 5 . 14 ) <!>, 


,^x{u) = exp(Y^ u^Xn (x). 
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In the sequel, when we will need to explieitly indieate the veetor fields defining <l> we will also use 
the notation: 

(5.15) ^e,x,x-{u) = and ^e,x,XH{u) its eomponents, 

and analogous notations will be used for the inverse map 

For every e > 0 and x,xq, in view of Theorem 13.91 the assoeiated ball box distanees reduee to: 

2 m u-\-m 2 u 

4(x,xo) = 

2=1 2=2m+l i=u+m-\-l 

For e = 0 and x,Xo we have 

n 

2=1 

5.4. Proof of the stability result. The sub-Laplaeian/heat operator assoeiated to this strueture is 

m+u 

Ka = dt-J2 
2=1 

where 

A = A®XI 

and I is the identity matrix of dimension u x u. We denote by F^^a the heat kernels of the eorre- 
sponding heat operators, and prove a lemma analogous to lemma 15^ for the lifted operator: 

Lemma 5.8. The fundamental solution of the operator L^ a> Is a kernel with local uniform 
exponential decay of order 2 with respect to e > 0 and A G according to definition d5.il) . 

Hence it belongs to the set S{2, d^, Moreover, e —)■ 0 one has 

(5.16) • --XldlT^^A ^ X,, • ■■X,,dlTA 

uniformly on compact sets, in a dominated way on all G. 

Proof The result for the limit operator Lo,a is well known and eontained for example in ifT^ . 
Henee we only have to estimate the fundamental solution of the operators a in terms of the one 
of To,A- In order to do so, we first define a ehange of variable on the Lie algebra: 

(5.17) T,{Xl) = X° for * = 1,..., + m 

Then from a fixed point 5 we apply the exponential map to induee on the Lie group a volume 
preserving change of variables. Using the notation introduced in (15.151) . we will denote 

G,-z :G^G, Ffix) = exp($-J_^^(^o)i(T)X°)(5) 

Since the distances are defined in terms of the exponential maps, this change of variables induces 
a relation between the distances do and dp. 

(5.18) dfix,Xo) = do{Ffix),Ffixo)). 








Analogously we also have 
(5.19) 
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Henee assertions (15.71) follow from the estimates of ro,A eontained for instance in ll54]| . Indeed 
the second inequality can be established as follows: 


T^^A{x,y,t) = To^A{Fe{x),F,{y),t) < Ca- 


df,{FA^),FAv)9 

3 Ca* 


= Ca- 


de{x,yY 

CAt 


Bo{F^{x),Vi)\ \B^{x,Vi)\ 


The proof of the first inequality in (15.71) and (15.81) is analogous, while (15.91) follows from the 
estimates of the fundamental solution contained in ( IfT^ l. 

The pointwise convergence (15.161) is also an immediate consequence of (15.181) and (15.191) . In 
order to prove the dominated convergence result we need to relate the distances do and On the 
other side, the change of variable (15.171) allows to express exponential coordinates u^, in terms of 
as follows: 


2m 


d.{x,xo) = E l“?l + E (l“? - luST-'"®)) 

SO that for all^ x, tq G (5 


2=1 


2=2m+l 


(5.20) 


do{x, Xo) -Co< de{x, Xo) < do{x, Xq) + Co 


where Co is independent of e. The latter and (15.81) imply that there is a constant Cs,k independent 
of e such that 


dn(x,yY 


1(5,-A-f ... A? r,A)(J,s. 91 < 


C\t 


\Bo{x,Vi)\ 

and this imply dominated convergence with respect to the e variable. 


□ 


In order to be able to conclude the proof of Proposition I5.21 we need to study the relation 
between the fundamental solutions ryi(x, y, t) and its lifting fo,A((a^, 0), {y, z),t), as well as the 
relation between re^(x, y, f)and f ^^^((x, 0), {y, z),t), 

Remark 5.9. We first note that for every / e x i?+) / can be identified with a C°° and 

bounded function defined on x and constant in the z— variables. Hence 


Consequently: 


f{.x,t) 


LaI = LaJ, L,,Af = Le,Af, 


r,,^((x,0 ,s), {y,z,t))dz 


Le,Afhj^s)dyds 


^This estimate indicates the well known fact that at large scale the Riemannian approximating distances are equiv¬ 
alent to the sub-Riemannian distance 

























30 


CAPOGNA AND CITTI 


From the definition of fundamental solution we can deduce that 
(5.21) 

VA{,x,y,t)= / Vo^A{,{.x,Q),{y,z),t)dz, and V^A{,x,y,t)= / V^^A{{.x,Q),{y,z),t)dz, 
Jg Jg 

for any x G G and f > 0. 


We conclude this section with the proof of the main result Proposition |5^ 


Proof. In view of the previous remanrk and (global) dominated convergence of the derivatives of 
T^^a to the corresponding derivatives of f o,a as e —)■ 0, we deduce that 



r,,A((x,0), {y,z),t)dz -)■ 



ro,A((a;,0), {y,z),t)dz 


as e —)■ 0. The Gaussian estimates of T^^a follow from the corresponding estimates on T^^a and the 
fact that in view of (15.201) . 


(5.22) 4((x, z), (xo, zo)) > do((x, z), (xq, zq)) - Cq > do(x, xq) + do(z, zq) - Cq > 


> de(x, Xo) + d^(z, Zo) - 3Co 

Indeed the latter shows that there exists a constant G > 0 depending only on G, ao such that for 
every x G G, 


de((x,^),(xo,zg)) 


The conclusion follows at once. 


d^(x,XQ) 


dz < Ce-t 


<ii(z,zo) 


df(x,XQ) 


t- dz < Ce ‘ 


>G 


□ 


5.5. Differential of the integral operator associated to T^. In this subsection we will show how 
to differentiate a functional F expressed as follows: 

F{f){x,t) = j T^^A{x,y,t)f{y,s)dyds. 

In order to do so, we will need to differentiate both with respect to x and to y, so that we 
will denote Xf^V^^^Aix, y, t) the derivative with respect to the variable x and Xf^V^^Aix, y, t) the 
derivative with respect to the variable y. 

Analogously, we will denote the derivative with the first variable of the lifted fundamental solu¬ 
tion 

Xf'^T^Aiix^w), {y,z),t). 

For e = 0, we will have by definition 

X°’"r,^((x, w), {y, z),t) = (X^ + Znr,A((F w), {y, z),t). 

The derivative with respect to the second variable will be denoted X^’^. If F is the Euclidean heat 
kernel, there is a simple relation between the derivative with respect to the two variables, indeed in 
this case T^^Aix, y, t) = T^^Aix — y, 0, t), so that 

(5.23) 2frr.,A(x, y, t) = -X,^’"r,,^(x, y, t). 
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Consequently for every funetion / e 

d^^F{f){x,t) = j V,^A{x,y,t)dyJ{y)dy. 

This is no more the ease in general Lie groups, or for Hormander veetor fields. However we will 
see that there is a relation between the two derivatives, whieh allows to prove the following: 

Proposition 5.10. Assume that f G C^(Hx]0, T[) in an open f2x]0, T[. For every x G iC CC 
H, for every i = 1, • • - m there exists the derivative XlF{f){x,t). Precisely there exist kernels 
R,,i(x,y,t) e £(2, such that 

X‘F{f) = 

/ m « 

^X"f^’*P,^i^h{x,y,t)f{y)dy- / R^^i{x,y,t)f{y)dy. 

h=i 

(Let us note explicitly that the term R^^i^h{x, y, t) plays the role of an error term). 

Proof We ean apply the lifting proeedure deseribed in seetions 5.2 and 5.3, and representing the 
fundamental solution as in (15.191) and (15.211) . we obtain the following expression for Fp. 

Fe{f) = J j f,,^((x,0), (|/,z),f)dz/(|/)di/= 

= J J ^o,AiFe{x,0),F^{y,z),t)dzf{y)dy. 

By differentiating with respeet to we get: 

(5.24) X^F,{f){x) = 11 - Zr)foAFe{x,0),F,{y,z),t)dzf{y)dy. 

Note that the family of veetors is independent of e and free of step r. Henee, by (see iTT^ . pag 
295, line 3 from below) for every i,j = 1, • • • m, sure exist families of indiees Rj, and polynomials 
pih homogeneous of degree > h sueh that: 

X°’^fo,A(T,i/,f) = 

m 

= Euf)* E xr’'(pii(e,(9))r„,^($,g.t))- 

i=l 

m 

E Xf')(p«(ej(!/)))ro,A(x,!/,(). 

j=l 

In partieular using this expression in the variable z alone, and integrating by parts we deduee 

11 ZrToAF.{x,0),FRy,z),t)dz = 0 
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We now call 

m 

RoA^^y^t) = Y 1 ^h^){pihi^Ay)))^o,Ai^^yR) 

J —1 h^Iij 

This kernel, being obtained by multiplication of ro,A(^, y, i) by a polynomial, has locally the same 
decay as ro,A(^, y, t)- In particular it is clear that the conditions l5.7ll5.8[|5.9l are satisfied uniformly 
with respect to e, since there is no dependence on e. As a consequence, if we set 

ReA^, y,t) = J RoAReA, 0), F^{y, z),t)dz 

then ReA^i y-i t) ^ ^^(2, d^, A Similarly we call 

Pe,i,A^^yR) = Y1 ^h^{piA^Ay))'^oA^^yR)) 

Now we use the fact that To^yi G S{2, d, together with the fact that pi^ is a polynomial of 

the degree equal of the order of to conclude that 

Ai,A^^ yR)^^ ( 2 , A mAvA 

It follows that, if we call 

Pe,i,A^^yR) = / h,i,h{RAA),Fe{y,z)R)d^ 


then Ax, y,t) G S{2,de,MAA) 

Plugging this expression into equation (15.241) we get 


(5.25) 


X'iF,{f)(x) = - Y.Xi''’rP..i.h(x,y,t)f{y)dy 


i=i 


ReA^^yR)fiy)dy- 


□ 


6. Stability of interior Schauder estimates 

In this section we will prove uniform estimates in spaces of Holder continuous functions and in 
Sobolev spaces for solutions of second order sub-elliptic differential equations in non divergence 
form 

n 

L^^au = dtu - ^ a\Ax, t)X-XA = 0, 

*j=i 

in a cylinder Q = x (0, T) that are stable as e —)■ 0. 

Indeed the proof of both estimates is largely based on the knowledge of a fundamental solution. 
Internal Schauder estimates for these type of operators are well known. We recall the results 
of Capogna and Han [?] for uniformly subelliptic operators, of Bramanti and Brandolini [fTTll for 
heat-type operators, and the results of Lunardi ll63 . and Polidoro and Di Francesco (SSI, and 
Gutierrez and Lanconelli ||4^ . which apply to a large class of squares of vector fields plus a drift 
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term. We also reeall [l64ll where uniform Sehauder estimates for a partieular elliptie approximation 
of subLaplaeians are proved. 

Here the novelty is due to the uniform eondition with respeet to e. This is aeeomplished by 
using the uniform Gaussian bounds established in in the previous seetion. This result extends to 
Hormander type operators the analogous assertion proved by Manfredini and the authors in |[T1 in 
the setting of Carnot Groups. 


6.1. Uniform Sehauder estimates. Let us start with the definition of elasses of Holder eontinuous 
funetions in this setting 


Definition 6.1. Let 0 < a < 1, Q C and u be defined on Q. We say that u e C^^xiQ) if 
there exists a positive eonstant M sueh that for every (x, f), (xq, to) ^ Q 

(6.1) \u{x,t) -u{xo,to)\ < M(i"((x,t), (xo,to)). 


We put 


\u\ 


c^xiQ) 


\u{x,t) - u{xo,to)\ 

= sup - 

d“((x,t), (xo,to)) 
k,a / 


sup \u\ 
Q 


Iterating this definition, if /c > 1 we say that u e C^'^xiQ) if for alH = 1 ,..., m X* e 


Where we have set = C^^xiQ)- 


The main results of this seetion, whieh generalizes to the Hormander veetor fields setting our 
previous result with Manfredini in [[Hi is 


Proposition 6.2. Let w be a smooth solution of = f on Q. Let K be a compact sets such 

that K C<Z Q, set 25 = do{K, dpQ) and denote by Ks the 5—tubular neighborhood of K. Assume 
that there exists a constant 67 > 0 such that 




e,X\ 


<C, 


for any e G (0,1). There exists a constant 67i > 0 depending on a, 67, 5, and the constants in 
Proposition \5.2\ but independent of e, such that 


\W\ \r^k+2, 

I ' Y 


\K) 






We will first eonsider to a eonstant eoeffioient operator, for whieh we will obtain a representation 
formula, then we will show how to obtain from this the elaimed result. 

Preeisely we will eonsider the eonstant eoeffioient/rozen operator: 


n 

Le,ixo,to) =dt-Yl aUxo,'to)X^X], 

Li=i 

where (xq, to) ^ Q- We explioitly note that for e > 0 fixed the operator is uniformly 

parabolio, so that its heat kernel oan be studied through standard singular integrals theory in the 
oorresponding Riemannian balls. 

As a direot oonsequenoe of the definition of fundamental solution one has the following repre¬ 
sentation formula 
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Lemma 6.3. Let w be a smooth solution to L^w = f in Q C For every cf) G C^{Q), 

(6.2) {w(p){x,t) 

= n^o4o)((^’ ^)) (w 0)(|/, T)dydT+ 

P n 

+ / r(^o,to)((a:, t), (y, r)) {jcj) + wL,(j) + 2 ^ T)X^wX](j)'^ {y, T)dydT, 

''Q i,j=i 

where we have denoted by the heat kernel for of L^^(^xoM)- 

Iterating the previous lemma we get the following 

Lemma 6.4. Let k E N and consider a k—tuple (ii ,... ,ik) G {1,.. . ,m}^. There exists a 
differential operator B of order k + 1, depending on horizontal derivatives of ajj of order at most 
k, such that 

n 

K ■■■K - A) = E (“« - “F^o, k))xi ■ ■ ■ XIX‘X‘ + B. 

*J = 1 

Proof. The proof ean be made by induetion. Indeed it is true with i? = 0 by definition if /c = 0: 


Le,(xo,to) — 'y ^ 

ij=l 

if it true for a fixed value of k then we have 




n 


B 


i,j=^ 


where 


B = io))X1 • ■ ■ X-,X-X- + A\^,B 


By the properties of B it follows that B is a differential operator of order k + 2, depending on 
horizontal derivatives of aof order at most k + 1. This eoneludes the proof. □ 


We can go back to our operator L and establish the following regularity results, differentiating 
twice the representation formula: 

Proposition 6.5. Let 0 < a < 1 and w be a smooth solution of L^w = / G CfxiQ) 
cylinder Q. Let K be a compact sets such that K CG Q, set 26 = do{K, dpQ) and denote by Ks 
the 5—tubular neighborhood of K. Assume that there exists a constant C > 0 such that for every 
e G (0,1) 

MjWci^iKs) < C. 

There exists a constant Ci > 0 depending on 6, a, C and the constants in Proposition 15.21 such 
that 
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Proof. The proof follows the outline of the standard ease, as in [|42l, and rests erueially on the 
Gaussian estimates proved in Proposition l5.2[ Choose a parabolie spherelli?^^^ <Z<Z K where <5 > 0 
will be fixed later and a eut-off funetion f G identieally 1 on i?e, 5/2 and eompaetly 

supported in B^ s- This implies that for some eonstant G > 0 depending only on G and do, 

|V, 0 | < C6-\ < C6-^, 

in Q. Now we represent the funetion wf through the formula [63] and take two derivatives in the 
direetion of the veetor fields. We remark onee more that the operator is uniformly elliptie due to 
the e—regularization, henee the differentiation under the integral ean ben eonsidered standard. As 
a eonsequenee for every multi-index I = (ii, * 2 ) ^ {1, • • •, and for every (xq, to) G one 
has: 

(6.3) XlX^^{wf){xo,to) 

(6.4) = / (|/,r))|(^o,4o) {L,,(xo,to) - L,) {w (j)){y,T)dydT+ 

Jq 

+ f Y‘,A7.rL,io)('. (!/. ^))(»,<.) ( f'l’ + “'i.'#' + 2 E 

•'e V i,,.i 

In order to study the Holder eontinuouity of the seeond derivatives, we note that the uniform Holder 
eontinuity of a-^-, and Proposition 15.21 ensure that the kemal satisfy the elassieal singular integral 
properties (see 0811): 


a‘jX‘wX‘4' I (j/, T)dydT. 


|ALYTL,)((^, i), (y,r)) - h), (y, t))| 


<le(xQ,y)^ 


< Cd'^{{x,t), {xo,to)) 


(r — to) 'T'AU-to) 

|B,( 0 , y/r - to)I 


with G > 0 independent of e. From here, proeeeding as in [|42l Theorem 2, Chapter 4], the first 
term in the right hand side of formula (16.31) ean be estimated as follows: 


(6.5) 


iyx)){Le - L,^(xo,to)){w (j)){y,T)dydT 




<Gi 




{Le Lf:(^xo,to)){'^f) 

= CiJ2\\{atj{xo,to) - atj{-))X]X]{w(j))\\^^^^^^^^^^ 
hj 

< Gi5" 11 alj 11 ) 11 11 ^.2,. , 


^That is a sphere in the group G = G x K in the pseudo-metric defined in (??). 
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where Ci, and Ci are stable as e —?■ 0. Similarly, if 0 is fixed, the Holder norm of the seeond term 
in the representation formula (16.31) is bounded by 


( 6 . 6 ) 


to), (y, T)){f(p{y, t) + wL())(p, r) + 2aljX^wX](f))dydT 




<C 2 

From (16.31) . (16.51) and (16.61) we deduee that 


67,, ,, 


\{Ks) 


\wc 


—672(5 I +672 


'Ct:xiBs) - " II llc^^xi^s) + ^2 11^1 lc);“(A5) 


67, 


Choosing 5 suffieiently small we prove the assertion on the fixed sphere The eonelusion fol¬ 
lows from a standard eovering argument. 

□ 


We ean now conclude the proof of proposition 16.21 

Proof. The proof is similar to the previous one for k = 1. We start by differentiating the represen¬ 
tation formula (16.21) along an arbitrary direction Xi^ 

(6-V) Xl{w(t)){x,f) 

( 6 . 8 ) = [ {y,T)) {L,^(xo,to) - Le) {w (l)){y,T)dydT+ 

Jq 

+ [ \f(l) + wL,(j) + 2^ 

V i,i=i 

Now we apply Theorem 15 . 101 and deduce that there exist kernels 

Pe,ii,h,{xQ,to){{x, f), {y, t)), Re,ii,{xo,to){{^i t), {y, 'P))) 

with the same decay of the fundamental solution such that 


afXlwXU {y, T)dydT. 


(6.9) 


Xl{w(t)){x,f) = 


/ m 

5 ^T’6,ii,h,(xo,to)((^)6), {y,T))Xl^y {Le,{xo,to) - Lf) (w (p){y,T)dydT- 

h=i 

- / Re,h,{xo,to){{x,t),{y,T)) - Le) {w f) {y, t) dy dx dy - 
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( 2 /, T)dydT. 


Using Lemma \6M this yields the existence of new kernels 
the behavior of a fundamental solution (and the same dependence on e) such that 


((x,t), (p,T))(^a^j - a^j(xo,to)^X^^^ ■ ■ ■ XlX^^XJ(wc/>)(p,T)dpdT 


n ' 


+ 


,(xq, to) 


{y,'r))B{w (t)){y,T)dpdT 


+ 


J e,hi,--- ,hfi,{xo,to) 

lower order terms 


n 

{{x, t), {p, t))XI ■ ■ ■ XI (^f(j){p, t) + wL,(t){p, t) + 2 aljX^wX](j)^dpdT+ 


i,j=^ 


where 0 is as in the proof of Proposition 16.51 and 5 is a differential operator of order k + 1. 
The conclusion follows by further differentiating the previous representation formula along two 
horizontal directions X^_^Xj^ and arguing as in the proof of Proposition [631 D 


7. Application I: Harnack inequalities eor degenerate parabolic quasilinear 

EQUATIONS HOLD UNIEORMLY IN e 

The results we have presented so far show that for any eo > 0, the 1—parameter family of met¬ 
ric spaces (Ai, d^) associated to the Riemannian approximations of a subRiemannian metric space 
{M,do), satisfy uniformly in e G [0, eo] the hypothesis in the definition of p-admissible structure in 
the sense of [j^ Theorem 13.1]. This class of metric measure spaces has a very rich analytic struc¬ 
ture (Sobolev-Poincare inequalities, John-Nirenberg lemma, ...) that allows for the development 
of a basic regularity theory for weak solutions of classes of nonlinear degenerate parabolic and 
elliptic PDE. In the following we will remind the reader of the definition and basic properties of 
p—admissible structures and sketch some of the main regularity results from the recent papers O 
and [[T9l . We will conclude the section with a sample application of these techniques to the global 
(in time) existence of solutions for a class of evolution equations that include the subRiemannian 
total variation flow [HI. 

7.1. Admissible ambient space geometry. Consider a smooth real manifold M and let p be 
a locally finite Borel measure on M which is absolutely continuous with respect the Lebesgue 
measure when represented in local charts. Let d{-, ■) : M x M ^ M+ denote the control distance 
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generated by a system of bounded, /i-measurable, Lipschitz veetor fields S = (Xi,..., Xm) on 
M. As in [in and [|45]| one needs to assume as basie hypothesis 

(7.1) the inelusion z : {M, chart) —)■ {M,d) is eontinuous, 

where we have denoted by (M, chart) the topology on M indueed by the Euelidean topology in 
via coordinate charts. For x & M and r > 0, set B{x,r) = {?/ G M : d{x,y) < r} and let 
\B{x, r)| denote the /i measure of B{x, r). In general, given a function u and a ball B = B{x, r) 
then Mb denotes the /i-average of u on the ball B = B{x,r). In view of (17.11) the closed metric ball 
.B is a compact set. 

Definition 7.1. Assume hypothesis (17.11) holds. Given 1 < p < cx), the triplet (M, /r, d) is said to 
define a p-admissible structure (in the sense of [|49l Theorem 13.1]) if for every compact subset K 
of M there exist constants Cd = CdC^, K), Cp = Cp{E, K) > 0, and R = i?(S, K) > 0, such 
that the following hold. 

(1) Doubling property: 

(D) |B(a:, 2r)| < CD\B{x,r)\ whenever a: G K and 0 < r < R. 

(2) Weak (1, p)-Poincare inequality: 

\u — upldp < Cpr I -r \E.uY'dp 

x,r) \J B{x,2r) 

whenever xeK,0<r<R,ue W^'^{B{x, 2r)). 

Theorems 11.11 and 11.21 yield the following 

Theorem 7.2. Let Xi,Xm be a family of Hormander vector fields in 12 C and denote by p 
Lebesgue measure. For each eo > 0 and e G [0, cq] denote by d^ the distance functions defined in 
Definition 12.61 For all e G [0, eo] andp > 1, the space (12,/r, dj is p—admissible, with constants 
Cp, and Cp independent of e. 

Other examples of p—admissible spaces are: 

• The classical setting: Xi = dp equals the n-dimensional Lebesgue measure, and 
E = {X,,...,Xm) = {d,„...,d,„). 

• Our setting is also sufficiently broad to include some non-smooth vector fields such as the 
Baouendi-Grushin frames, e.g., consider, for 7 > 1 and {x,y) G M^, the vector fields 
Xi = dx and X 2 = \x\'^dy. Unless 7 is a positive even integer these vector fields fail to 
satisfy Hormander’s finite rank hypothesis. However, the doubling inequality as well as the 
Poincare inequality hold and have been used in the work of Franchi and Lanconelli |[39ll to 
establish Harnack inequalities for linear equations. 

• Consider a smooth manifold Xi endowed with a complete Riemannian metric g. Let p 
denote the Riemann volume measure, and by S denote a p—orthonormal frame. If the Ricci 
curvature is bounded from below (Ricci > —Kg) then one has a 2—admissible structure. 
In fact, in this setting the Poincare inequality follows from Buser’s inequality while the 
doubling condition is a consequence of the Bishop-Gromov comparison principle. If K = 
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0, i.e. the Ricci tensor is non-negative, then these assumptions holds globally and so does 
the Harnack inequality. 

Spaces with a p-admissible structure support a homogenous structure in the sense of Coifman and 
Weiss If29]| . 

Lemma 7.3. Let (M, ju, d) be a p-admissible structure for some p > 1, LL a bounded open set in 
M. and set K = fl. If x E K and 0 < s < r < R, then the following holds. 

(1) There exists a constant N = N{Cd) > 0 , called homogeneous dimension of K with 
respect to (S, d, p), such that \B{x, r)| < Cdt~^\B{x, Tr)\, for all 0 < r < 1. 

(2) There exists a continuous function 0 G Cq{B{x, r)) n W^’^{B{x, r)) and a constant C = 
C(S, K) > 0, such that f = lin B{x, s) and |S0| < C/{r — s), 0 < f < 1. 

(3) Metric balls have the so called 5—annular decay property, i.e., there exists S = S^Cd) G 
( 0 , 1 ], such that 

\B{x, r) \ B{x, (1 — e)r)| < Ce^\B{x, r)|, 
whenever 0 < e < 1. 


Proof. Statement (1) follows from dD]) by a standard iteration argument. Statement (2) is proved 
in [|^ Theorem 1.5]. Statement (3) follows from [[T3l Corollary 2.2], since we have a Carnot- 
Caratheodory space. Furthermore, 5 depends only owCd. □ 


Given Vt <Z M, open, and 1 < p < oo, we let = {u E p) : XiU E p),i = 

1,..., m} denote the horizontal Sobolev space, and we let C be the closure0of the space 
of functions with compact (distributional) support in the norm ||m||i p = ||m||p + ||Sm||p with 
respect to p. In the following we will omit p in the notation for Lebesgue and Sobolev spaces. 
Given R < t 2 , and 1 < p < oo, we let x (^ 1 ,^ 2 ) and we let LP(ti,t 2 ; ti < t 2 , 

denote the parabolic Sobolev space of real-valued functions defined on such that for almost 
every t,ti<t< t 2 , the function x —)■ u{x, t) belongs to lF’i’^(f 2 ) and 


\u\ 






('*2 




1/p 


' ti J Q 


{\u{x,t)f + \Eu{x,t)f)dpdt] < 00 . 


The spaces LP{ti,t 2 ', IFi’Q(f2)) is defined analogously. We let ^ 2 ; consist of real¬ 

valued functions p G L^{ti,t2] L^iVt)) such that the weak derivative dtp{x, t) exists and belongs to 
L^(ti,t 2 ] LP{fl)). Consider the set of functions (f>, (f> E t 2 ] L^(f2)), such that the functions 


f —> f \(j){x,t)fdp{x) and t ^ f \dt(j){x,t)fdp{x), 

Jn Jn 

have compact support in (ti, ^ 2 )- We let WQ’^{ti,t2] L^{Q)) denote the closure of this space under 
the norm in ^ 2 ; Lp{Q)). 

From ll49l Corollary 9.5] one can see that the metric balls B{xo, r) are John domains. Conse¬ 
quently, dni), (0, and [|43 Theorem 9.7] yield Sobolev-Poincare inequality. 


^Here we avoid the issue “iJ = W". This is studied in detail in g3l . ED, mi, iOl, MB and (EH. 
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Lemma 7.4. Let B{xo,r) gQ, 0<r<R, l<p<oo. 
C{Cd,Cp,p) > 1 such that for every u G W^^{B{xQ,r)), 

/ r \ r 


\u-UBf^dp\ <CrP 


There exists a constant C = 




\J B{xo,r) / J B{xo,r) 

where ub denotes the p average ofu over B{xo, r), and where 1 < n < N/ {N — p), ifl <p < N, 
and 1 < n < oo, ifp > N. Moreover, 


1/K 

\ufPdp ] < CrP 


B{xo,r) 


/ 


J B{xQ,r) 


lufdp, 


whenever u G W^’Q{B{xQ,r)). 


7.2. Quasilinear degenerate parabolic PDE. In this section we list some recent results concern¬ 
ing regularity of weak solutions of certain nonlinear, degenerate parabolic PDE in spaces (M, p, d) 
that are p-admissible for some p G [2, oo). If p = 2 we can allow lower order terms, but at present 
this is not yet established for p > 2. Given a domain (i.e., an open, connected set) D C M, and 
T > 0 we set Up = 12 x (0, T). For a function u : 12^ —?■ M, and 1 < p,qwe define the norms 


(7.2) 




\ufdx)i‘dt 


and the corresponding Lebesgue spaces L^’^iVtT) = -^^([0, T], Lp( 12)). We will say that A, B are 
admissible symbols (in Lip) if the following holds: 


(i) (x, t) —)■ A{x, t, u, ^), B{x, t,uA) are measurable for every (m, 0 £ ® x M™, 

(ii) (m, 0 A{x, t,uA), B{x, t,uA) are continuous for almost every (x, t) G Lip, 

(iii) • For p = 2: There exist constants a, a > 0 and functions b,c,e,f,h G LP’‘^{Q) with 

p > 2, and q given by ^ ^ and functions d, p G L°‘’^{Q) with 1 < a and /) given 

by ^ i < 1 such that for a.e. (x, t) G 12^ and ^ G one has 


(7.3) 


(7.4) 


^A^{xA,uA)^i 

i=l 

\A{x,t,uA)\ 

\B{x,t,uA)\ 


>a|ep-6V-/^ 

< d|.^| + e\u\ + h, 

< c\^\ +d\u\ +g. 


In view of the conditions on p,q, a, (3 there exists 6 > 0 such that 

2 N 1 1-9 

^ - 1-9 2p q - 2 

1 . 2V 1 

a > -- and- \- — < 1 — 9. 

-1-9 2a /3 - 
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We say that a constant depends on the structure conditions (17.31) . if it depends only on0 

||6||, ||c||, ||d||, ||e||, ll/ll, ll^ll, ||/i||,iV,0, 

and is uniformly bounded if these quantities are so. 

• For p > 2 we will only consider S = 0 and ask that the following bounds 

(7.5) A{x,t,u,^) ■ ^ > Ao\^\P, \A{x,t,uA)\ < Ai\^\P~\ 

hold for every (uA) G M x and almost every (x, t) G fir- 

Aq and .4,1 are called the structural constants of A. If A and A are both admissible symbols, with 
the same structural constants Aq and Ai, then we say that the symbols are structurally similar. 

Let if be a domain in Ai x M. We say that the function u : if —)■ M is a weak solution to 


m 

(7.6) dtu{x, t) = La,pU = — X*Ai{x, t, u, Eu) + B{x, t, u, Eu), 

i=l 


in if, where X* is the formal adjoint w.r.t. dp, if whenever ^ E for some domain fl C Ai, 
u e LP{ti,t 2 ', lA’s’^(^)) 


(7.7) [ u^dpdt 

J U J Q 


f*t2 


dt 


A{x, t, u, Eu) ■ Er] dpdt + 


f*t2 


B{x, t, u, Eu)r} dpdt = 0, 


for every test function 

A function u is a weak super-solution (sub-solution) to (17.61) in E if whenever E for some 

domain fl C Ai, we have u G L^(ii, ^ 2 ; kF^’^(n)), and the left hand side of (17.71) is non-negative 
(non-positive) for all non-negative test functions kFQ’^(ti, ^ 2 ; Fl L^(ii, ^ 2 ; o(^))- 

The main results in O and [|T^ can be summarized in the following theorem. 


Theorem 7.5. Let (Ai, p, d) be a p-admissible structure for some fixed p G [2, 00 ). For a bounded 
open subset fl C Ai, let u be a non-negative, weak solution to (17.61) in an open set containing the 
cylinder fl x [0, Tq] and assume that the structure conditions (17.51) are satisfied. 

• For p = 2 and for any subcylinder Q^p = B{x,3p) x (f — 9p‘^,t) C Q there exists a 
constant C > 9 depending on Cd, Cl, Cp, the structure conditions (17.31) and on p such 
that 

(7.8) supu < (7inf(M + p®/c), 

Q- Q+ 

where 

(7.9) = B{x,p) X {t — p^,t) and Q~ = B{x, p) x {t — 8p‘^,i — 7p^) 

9 > 0 is defined as in (17.41) . and we have let k = 11/| | + | |p| | + | |/r| |- 


^ The 11 • 11 norms are in the appropriate or classes 
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• For p > 2: Assuming B = 0, there exist constants Ci, C 2 , C 3 > 1, depending only on S, 
CoiCp, Aq, Ai^p, such that for almost all {xqAo) G x [0,To], the following holds: If 
m(xo, to) > 0, and i/O < r < -R(S, Q) (from Definition \7.1\) is sufficiently small so that 


B{xo, 8 r) C and (to - C'iM(a;o, to)^ to + Cim(xo, to)^ C (0,To), 
then 


u{xo,to) < C 2 infM, 
Q 


where 


Q = B(xo,r) X (^o + ^C3u(xo,tof + C3u(xo,toY 


Furthermore, the constants Ci, C 2 , C 3 can be chosen independently ofp as p ^ 2. 

We conclude this section with a corollary of the proof in lfT9li rLemma 3.6], a weak Harnack 
inequality that plays an important role in the proof of the regularity of the mean curvature flow 
for graphs over certain Lie groups established in ifT^ . Consider a weak supersolution w G 
U‘(tiAA of the linear equation 


(7.10) - dtw - '^X*(aij(xA)Xjw) = g{xA), 

i=\ 


with ti, t2, ft as defined above. Assume the coercivity hypothesis 


(7.11) 


m 

d{i)=l *,i=l d{i)=l 


for a.e. (x, t) and all f G M™, for a suitable constant A. 

Proposition 7.6. Let {Ai,p,d) be a 2-admissible structure. For a bounded open subset Ft C 
Xi, let u be a non-negative, weak supersolution to (17.101) in an open set containing the cylinder 
Ft X [0, To] and assume that conditions (17.1 II) are satisfied. For any subcylinder = B(x, 3p) x 
(t — 9p^, t) G Q there exists a constant (C > 0 depending on Cp, Cp, Cp, the structure conditions 
(17.31) and on p such that 

(7.12) f w dxdt < C(ini w + sup \g\p‘^), 

\Q \ JQ- < 5 + Q+ 


with Q~^, Q as defined in (17.91) . 
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8. Application II: Regularity for quasilinear sub elliptic PDE through 

Riemannian approximation 


As an illustration of the usefulness of the uniform estimates established above, in this seetion 
we want to briefly sketeh the strategy used in [fT^ and [fTTll . where the Riemannian approximation 
seheme is used to establish regularity for the graph solutions of the Total Variation flow 


( 8 . 1 ) 


du _ ^ / XiU 

^ Vv^TTTv^ 


and for the graphieal solutions of the mean eurvature flow 


( 8 . 2 ) 


du 

'dt 


i=l 


XjU 

Vi + IVomP 


In both oases D C G is a bounded open set, with G is a Lie group, free up to step two, but not 
necessarily nilpotent. 

We will consider solutions arising as limits of solutions of the analogue Riemannian flows, i.e. 


(8.3) 

and 

(8.4) 


du^ 

dt 


Xfu, 


n 

i=l 


for T G D, t> 0, 


dUf^ 

dt 


wx 



for T G D, t > 0, 


where, is the mean curvature of the graph of Me(-, f) and 

n . . 

(8.5) W^ = l + |V.t<.P = 1 + and 

for all G M”. 

The main results in [[T^ and [fTTll concern long time existence of solutions of the initial value 
problems 

(8 6) / ^ in Q = D X (0, T) f dtu ^ = in Q = D x (0, T) 

^ I Me = </? on dpQ, I Me = on dpQ, 

with dpQ = {Q X {t = 0}) U {dQ x (0, T)) denoting the parabolic boundary of Q. 


Theorem 8.1. Let G be a Lie group of step two, Vt <Z G a bounded, open, convex set (in a sense 
to be defined later) and ip G G^(r2). There exists unique solutions Me G G°°(D x (0, cx))) fl 
L°°((0, CX)), G^(r2)) of the two initial value problems in (18.61) . and for each k E N and K CC Q, 
there exists G^ = Gfe(G, p, k, K,Vl) > 0 not depending on e such that 

(8.7) I |We| |c'=(A') ^ C'fc. 
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Corollary 8.2. Under the assumptions of Theorem l&il as e ^ ft the solutions of either flow 
converge uniformly (with all theirs derivatives) on compact subsets of Q to the unique, smooth 
solution of the corresponding sub-Riemannian flow in Vt x (0, cxo) with initial data p. 


The proof of this result rests erucially on the estimates established in this paper. In the following 
we list the main steps. First of all we note that in view of the short time existence result in the 
Riemannian setting we can assume that locally are smooth both in time and space. 

(1) Interior gradient bounds. Denote by right XI the left invariant frame corresponding 
to X^s and observe that these two frames commute. For both flows, consider solutions 
Me e C^iQ) and denote vq = dtu^, Vi = X^u^ for i = i, ..., n. Then for every /i = 0,..., n 
one has that Vh is a solution of 


( 8 . 8 ) 


where 



1 

TTTfF 


for the total variation flow, while 



i + \^\y 


4(0 = 4 


i + iep’ 


for the mean curvature flow. The weak parabolic maximum principle yields that there exists 
C = C{G, 119?| |c 2 (n)) > 0 such that for every compact subset K CC Tl one has 


sup |VlMe| < SUp(|VlMe| + \dtUe\), 

Kx[0,T) dpQ 


where Vi is the full ^fi—Riemannian gradient. This yields the desired unform interior 
gradient bounds. This argument works in any Lie group, with no restrictions on the step. 

(2) Global gradient bounds. The proof of the boundary gradient estimates is more deli¬ 
cate and depends crucially on the geometry of the space. In particular the argument we 
outline here only holds in step two groups G and for domains C G that are locally Eu¬ 
clidean convex when expressed in the Rothschild-Stein preferred coordinates introduced in 
(15.1 II) . In ifTSll we use the Rothschild-Stein osculation Theorem 15.71 to construct a rather 
explicit barrier function at any boundary point and then to conclude we apply the compar¬ 
ison principle lfT4l Theorem 3.3]. This argument also shows that the solutions Vh to (18.81) 
are bounded. 

(3) Harnack inequalities and G^’°‘ estimates. We have proved in Theorem [LTI and Theorem 
11.21 that (G, dfj is a 2—admissible geometry in the sense of Definition 17.II with Doubling 
and Poincare constants uniform in e > 0. As a consequence we can apply the Harnack 
inequalities in Theorem 17.51 and Proposition 17.61 to the bounded solutions Vh of (18.81) . thus 
yielding the G^’°‘ uniform interior estimates. 

(4) Schauder estimates and bigber order estimates The uniform Gaussian estimates and 
Schauder estimates in Theorem 11.41 applied to (18.81) yield the higher order estimates and 
conclude the proof. 
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